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ABSTRACT 


The  movement  and  dispersion  of  vehicles  (cars,  trucks 
and  buses)  as  they  travel  between  traffic  signals  is  an 
important  practical  consideration  in  both  the  design  and 
analysis  of  fixed-time  signal  co-ordination  timing  plans. 
Knowing  the  arrival  pattern  of  a  platoon  of  vehicles 
released  from  a  previous  set  of  traffic  signals  on  a  road  is 
a  requirement  when  calculating  the  signal  offset  and  setting 
of  a  subsequent  set  of  traffic  signals  on  the  same  road. 
Platoon  dispersion  models  provide  some  insight  into  this 
complex  phenomenon. 

Among  the  various  studies  in  platoon  dispersion 
(Lighthill's  and  Witham's  wave  theory,  Pacey's  diffusion 
theory  and  Robertson's  recurrence  model),  the  research  by 
Robertson  represents  a  major  advancement.  This  is  indicated 
by  the  incorporation  of  his  platoon  dispersion  model  in  such 
signal  optimization  programs  as  TRANSYT  (TRAffic  Network 
StudY  Tool)  and  the  GLC  (Greater  London  Council)  Combination 
Method.  Robertson's  platoon  dispersion  model  uses  a 
recurrence  exponential  smoothing  equation  to  predict  a 
downstream  arrival  pattern  (downstream  cyclic  flow  profile) 
from  an  upstream  arrival  pattern  (upstream  cyclic  flow 
profile).  The  smoothing  factor  is  a  function  of  the  average 
travel  time  between  the  upstream  and  downstream  cyclic  flow 
profiles,  and  two  empirical  parameter:  alpha  and  beta. 
Robertson  suggests  average  numerical  values  for  these 
parameters,  but  has  qualified  these  values  by  stating  that 
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the  parameters  may  vary  with  the  site  characteristics. 


This  thesis  concentrates  on  calibrating  Robertson's 
Platoon  Dispersion  Model  from  the  persective  of  vehicle 
delays,  number  of  stops,  and  performance  index  for  traffic 
platoons  in  Edmonton,  Alberta.  Actual  real-time  summer  and 
winter  platoon  data  was  collected  with  a  Micro-computer 
Cyclic  Flow  Profile  Monitor. 

The  analysis  indicated  that  the  platoon  dispersion 
predicted  by  Robertson's  model  using  his  recommended 
empirical  parameters  could  be  in  error  when  compared  with 
the  observed  data.  The  impact  of  this  error  on  signal 
performance  as  quantified  by  delay,  stops  and  performance 
index  was  significant.  The  analysis  revealed  that  improper 
empirical  parameters  and  the  specific  form  of  Robertson's 
model  contributed  to  the  error.  Robertson's  model  requires 
the  flow  value  in  the  previous  downstream  interval  to 
predict  the  value  in  the  present  interval.  This  information 
is  not  available  in  computing  the  flow  during  the  first 
iteration  and  is  assumed  to  take  a  zero  value.  Since  traffic 
platoons  between  fixed-time  signals  are  cyclic  in  nature, 
and  if  steady  state  conditions  exist,  the  flow  for  an 
interval  will  be  the  same  over  all  cycle  lengths.  Based  on 
this  information  an  initialization  equation  was  derived.  For 
best  results,  the  research  indicates  that  calibration  of  the 
individual  values  of  alpha  and  beta  along  with  a  modified 
form  of  Robertson's  model  is  essential. 


Based  on  the  survey  results  there  is  some  evidence  to 
suggest  that  the  individual  values  of  alpha  and  beta  are 
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1.  INTRODUCTION 


1 . 1  Problem  Statement 

The  movement  and  dispersion  of  vehicles  (cars,  buses 
and  trucks)  as  they  travel  between  traffic  signals  is  an 
important  practical  consideration  in  both  the  design  and 
analysis  of  fixed-time  signal  co-ordination  timing  plans. 
Knowing  the  arrival  pattern  of  a  platoon  (an  analysis  of 
vehicular  flow  rates  with  a  period  equal  to  the  cycle 
length,  which  indicates  a  higher  rate  of  random  flow,  or 
which  some  type  of  regularity  or  pattern  occurs,  due  to  an 
upstream  signal,  may  be  defined  as  a  platoon)  of  vehicles 
released  from  a  previous  set  of  traffic  signals  on  a  road  is 
a  requirement  when  calculating  the  signal  offset  and 
settings  of  a  subsequent  set  of  traffic  signals  on  the  same 
road.  The  application  of  this  relationship  is  illustrated  in 
Figure  1.1  and  described  below. 

If  the  cumulative  vehicular  demand  function  ( i . e . , 
number  of  vehicles  arriving  at  the  signal  approach)  and  the 
cumulative  vehicular  service  function  ( i  .  e .  ,  number  of 
vehicles  discharging  at  the  signal  stop  line)  are  plotted 
against  time  on  the  same  graph,  the  following  performance 
indicators  can  be  quantified  : 

1.  The  area  between  the  two  cumulative  functions  represents 
the  total  vehicular  delay, 

2.  A  vertical  line  between  the  two  cumulative  functions 
represents  the  vehicular  queue  or  the  number  of  vehicles 
stopped  at  a  given  time,  and 
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FIGURE  1.1 


REPRESENTATION  OF  QUEUING  AT  A  SIGNALIZED  INTERSECTION 
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3.  A  horizontal  line  between  the  two  cumulative  functions 
represents  the  vehicular  delay  for  a  given  vehicle. 

The  effect  of  varying  the  offset  time  between  the  two 
co-ordinated  traffic  signals  at  the  end  of  the  roadway  link 
will  be  to  move  the  cumulative  service  function  relative  to 
the  cumulative  demand  function;  thus  the  fundamental 
relationship  between  delay/stops/performance  index  and 
signal  offset  can  be  obtained. 

Central  to  this  fundamental  concept  is  the  observation 
that  the  cumulative  demand  function  is  influenced  by  the 
manner  in  which  a  platoon  of  vehicles  moves  along  a  roadway 
link.  Vehicles  in  the  platoon  will  attain  different  speeds 
due  to  a  variety  of  factors:  vehicle  characteristics,  driver 
behaviour,  road  surface  characteristics,  interference  from 
pedestrians,  roadside  activities  and  side  street  traffic. 
Ignoring  the  effect  of  this  phenomenon  would  reduce  the 
effectiveness  of  signal  coordination  design.  Platoon 
dispersion  models  provide  some  insight  into  this  complex 
phenomenon . 

Among  the  various  studies  in  platoon  dispersion,  the 
research  by  Lighthill  and  Whi tham ( 1 , 2 ) ,  Pacey(3)  and 
Robert son ( 1 2 )  represent  major  advancements  in  the 
development  of  platoon  dispersion  models.  The  work  by 
Robertson  is  by  far  the  most  significant  in  terms  of 
practical  application,  as  indicated  by  the  incorporation  of 
his  platoon  dispersion  model  in  such  signal  optimization 
programs  as  TRANSYT(13)  and  the  GLC  Combination  Method(23). 
The  effectiveness  of  these  programs  for  both  conventional 


. 


1:' 


4 


and  computerized  signal  control  systems  has  been 
demonstrated  in  a  number  of  field  appl icat i ons ( 2 4 , 2 E ) .  The 
significance  of  Robertson’s  research  is  partially  due  to  his 
intoduction  of  the  Cyclic  Flow  Profile  Concept!2*). 
Robertson’s  platoon  dispersion  model  is  based  upon  the 
theory  that  the  platoon  of  vehicles  leaving  an  upstream 
intersection  is  formed  by  the  queue  of  vehicles  waiting  at 
the  start  of  green,  departing  at  the  saturation  flow  rate 
and  dispersing  as  the  platoon  travels  along  the  segment. 
This  concept  satisfies  a  practical  need  in  that  it 
represents  a  compromise  between  using  micro-scopic  models 
(simulation  of  individual  vehicular  behaviour)  and 
macro-scopic  models  (relationships  between  vehicle  flow 
rates  and  speed) . 

Continued  research  in  the  area  of  platoon  dispersion 
has  shown  that  the  Robertson's  model  can  provide  a 
reasonable  fit  with  observed  data,  but  the  accuracy  of  the 
model  can  be  enhanced  significantly  by  calibrating  to  local 
conditions  and/or  improving  on  the  structure  of  the  model. 
Table  1.1  provides  a  comprehensive  and  chronological  summary 
of  published  calibrated  empirical  parameters  in  Robertson’s 
platoon  dispersion  model.  As  shown  in  Table  1.1,  there  are  a 
significant  range  of  empirical  parameters  which  result  in  a 
good  prediction  of  the  dispersion  characteristics  of 
observed  traffic  platoons.  Lam  (22)  indicated  that  a  lack  of 
calibration  may  result  in  the  prediction  of  minimum  delay 
estimates  which  are  in  error  in  the  range  of  5  -  20%.  It  is 
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TABLE  1 . 1 


COMPREHENSIVE  CHRONOLOGICAL  SUMMARY  OF  PUBLISHED  CALIBRATED 
EMPIRICAL  PARAMETERS  IN  ROBERTSON T  S  PLATOON  DISPERSION  MODEL 
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this  variablility  in  the  empirical  parameters  and  the  impact 
of  these  empirical  paramaters  on  such  signal  performance 
indicators  as  delay,  stops  and  performance  index  that  has 
prompted  this  research  for  the  Edmonton,  Alberta 
environment.  More  specifically,  this  research  effort 
concentrates  on  the  dispersion  of  traffic  along  roadway 
links  typically  in  the  Central  Business  District. 

In  addition,  cyclic  flow  profiles  provide  the 
transportation  engineer  with  a  great  deal  of  additonal 
information  relevant  to  signal  design  and  analaysis.  Such 
profiles  allow:  review  of  average  flow  conditions, 
identification  of  major/minor  road  inflows,  the  benefits  of 
signal  co-ordination,  estimates  of  saturation  flow, 
estimates  of  spare  capacity,  estimates  of  start  and  end 
lags,  estimates  of  platoon  speed  and  amount  of  dispersion, 
incident  detection  poss ibl i t i es ,  verification  of  signal 
settings  by  TRANSYT,  possibilities  for  an  international 
exchange  standard,  determination  of  preferential  street 
treatment  and  ranking  of  the  importance  of  network  links. 
They  also  point  to  the  need  for  coordination,  and  aid  in 
computer  network  condensation.  The  use  and  understanding  of 
this  additional  information  is  essential  and  a  pre-requisite 
to  good  signal  design  and  analysis. 


8 


1.2  Thesis  Objective 

The  primary  objective  of  this  thesis  is  to  calibrate 
Robertson's  Platoon  Dispersion  Model  for  data  in  Edmonton, 
Alberta.  Data  was  based  on  real-time  traffic  platoons  during 
summer  and  winter  conditions.  This  data  base  was  in  response 
to  the  practical  problems  associated  with  winter  driving 
conditions  in  Alberta  (33)  and  more  generally  across  Canada. 
To  achieve  the  objective  two  tasks  were  undertaken.  The 
first  task  involved  the  validation  of  the  basic  structure  of 
the  model  as  proposed  by  Robertson  and/or  some  modification 
to  its  basic  form  based  on  a  mathematical  extension.  The 
second  task  involved  the  calibration  of  the  model(s)  to 
several  "typical"  roadway  links  which  characterize  arterials 
in  the  Central  Business  District.  The  research  objective 
requires  the  following  outputs: 

1.  Documentation  of  the  difference  in  signal  performance 
characteristics  such  as  delays,  stops  and  performance 
index  between  a  calibrated  and  not  calibrated  platoon 
dispersion  model. 


2.  Presentation  of  design  guidelines  which  will  assist  in 
the  selection  of  the  most  appropriate  empirical 
parameters . 

A  secondary  objective  of  this  reseach  is  to  translate 
the  experience  obtained  during  the  course  of  this  reseach 
into  a  series  of  practical  applications  as  outlined  in 


Sec t i on  1.1 
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1.3  Scope  Of  The  Research 

In  organizing  this  research  project  a  number  of 
activities  are  defined  which  will  dictate  the  course  of 
action  in  meeting  the  research  objectives.  The  following 
considerations  will  assist  in  limiting  the  research  to  a 
manageable  level: 

1.  This  study  will  not  attempt  to  verify  and  calibrate  any 
platoon  dispersion  model  other  than  that  suggested  by 
Robertson . 


2.  All  street  segments  to  be  surveyed  will  be  operated 
under  a  fixed  time  signal  control  strategy. 


3.  The  selected  survey  locations  (total  of  6  measurements 
during  the  A.M.  peak  hour)  will  represent  typical 
arterial  roadways  within  or  leading  to  the  Central 
Business  District  (CBD)  in  Edmonton,  Alberta,  Canada. 


4.  Dispersion  analysis  will  be  restricted  to  link  lengths 
which  are  typical  of  the  CBD  environment. 


1.4  Organization  Of  The  Thesis  Document 

The  introduction  along  with  the  problem  statement, 
thesis  objective,  scope  of  the  research  and  organization  of 
the  thesis  document  are  given  in  Chapter  1. 

Chapter  2  reviews,  in  detail,  past  research  efforts  in 
the  area  of  platoon  dispersion.  Thoeretical  aspects  of 
Lighthill’s  and  Whitham's  Wave  Theory,  Pacey's  Platoon 
Diffusion  Theory  and  Robertson's  Recurrence  Model  are 
presented.  Mathematical  considerations  in  Robertson's 
Recurrence  Model  are  critically  discussed. 
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Chapter 

3  discusses  the  data 

collection 

efforts  with 

respect  to 

winter  and  summer  platoon 

data  and  speed/travel 

time 

data.  I 

n  addition,  the  survey 

site 

s  and 

their  observed 

A.M. 

peak 

hour  characteristics 

and 

site 

conditions  are 

presented . 

Chapter  4  presents  the  basic  analytical  tools  that  were 
used  in  the  analysis  process.  More  specifically,  three 
computing  programs  were  employed  to  aid  in  the  determination 
of  the  "best"  empirical  parameters. 

Chapter  5  presents  a  detailed  summation  of  all  analysis 
undertaken.  The  results  of  the  calibrated  process  are 
discussed  from  the  perspective  of  delay,  stops  and 
performance  index.  Design  guidelines  are  also  presented  to 
assist  in  the  selection  of  the  most  appropriate  empirical 
parameters  in  Robertson's  Model.  In  addition,  practical 
applications  of  observed  traffic  platoons  are  presented. 

Chapter  6  presents  the  conclusions,  recommendations  and 
further  research  suggestions  arising  from  this  research 


effort 


2.  PLATOON  DISPERSION  MODELS 


2.1  Overview  Of  Past  Research  Efforts 

The  amount  of  theoretical  and  applied  research  into  the 
behaviour  of  traffic  platoons  and  effects  of  their 
dispersion  on  traffic  control  has  increased  steadily  since 
1  955. 


One  of  the  earliest  (1955)  research  works  related  to 
platoon  dispersion  was  carried  out  by  Lighthill  and 
Whitham( 1 , 2 ) .  These  scientists  presented  their  kinematic 
wave  theory  in  two  papers,  the  former  giving  a  more 
mathematical  treatment  and  dealing  specifically  with  flood 
movement  in  long  rivers,  the  latter  giving  a  more 
descriptive  treatment  and  dealing  with  the  flow  of  traffic 
on  long  crowded  roads.  In  the  two  papers,  attention  was 
drawn  to  waves  which  can  be  expressed  in  a  one-dimensional 
flow  system.  These  waves  exist  if,  to  a  sufficient 
approximation,  there  is  a  functional  relationship  between 
the  flow  or  volume,  the  concentration  or  density  and  the 
distance  along  the  roadway  link.  On  the  assumption  of  a 
functional  relationship,  the  wave  theory  follows  from  the 
equation  of  continuity  alone.  Accordingly,  the  waves  are 
described  as  "kinematic",  as  opposed  to  "dynamic"  waves 
which  depend  on  Newton's  second  law  of  motion  and  some  other 
assumptions  relating  stress  to  displacement,  strain,  or 
curvature.  Kinematic  waves  do  not  disperse  as  other  waves 
do,  but  they  suffer  a  change  in  form,  due  to  the  dependence 
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of  the  wave  speed  on  the  flow  carried  by  the  waves. 
Accordingly,  kinematic  wave  forms  may  develop  di scont inut ies 
due  to  the  overtaking  of  slower  waves  by  faster  ones.  These 
were  described  as  "shock  waves"  since  they  were  formed  in 
the  same  way  as  shock  waves  in  a  gas.  From  a  plot  of  both 
the  continuous  kinematic  waves  and  shock  waves,  the  cyclic 
flow  profile,  and  thus  the  platoon  structure  at  any  point 
along  the  roadway  section,  could  be  derived. 

Pacey(3),  in  1956,  attempted  to  use  Lighthill  and 
Whitham's  kinematic  wave  theory  to  predict  platoon 
dispersion,  and  also  developed  his  own  platoon  prediction 
model  by  treating  platoon  dispersion  as  a  fluid  diffusion 
process.  The  fluid  diffusion  theory  assumes  that  the  only 
changes  in  shape  of  the  bunch  of  traffic  released  from  the 
signals  arise  from  differences  in  speed  between  vehicles  in 
the  bunch;  there  is  no  interference  with  overtaking,  and  any 
given  vehicle  in  a  platoon  proceeds  with  a  constant  velocity 
(normally  distributed  about  a  mean  speed)  irrespective  of 
the  number  or  distribution  of  the  vehicles  on  the  road. 

Lewis(4),  in  1958,  carried  out  a  number  of  surveys 
downstream  from  a  signal  on  a  semi-expressway  facility  and 
found  that  signal  coordination  is  feasible  for  a  link 
distance  of  up  to  1.0  kilometre  (1,000  metres)  or  0.65  miles 
(3,430  feet) . 

Gerlough(5),  in  1961,  suggested  that  an  analogy  exists 
between  traffic  dynamics  and  wave  mechanics  which  could  be 
used  to  describe  the  dissipation  of  traffic  platoons  when 
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released  from  a  traffic  signal. 

Graham  and  Chenu(6),  in  1962,  studied  the  phenomenon  of 
platoon  dispersion  at  various  points  downstream  from  a 
signal.  Field  studies  confirmed  that  at  points  405,  805, 
1210,  and  1610  metres  (1/4,  1/2,  3/4  and  1.0  mile),  the 
percent  of  vehicles  remaining  in  the  platoon  was  91,  85,  80 
and  77%,  respectively.  This  applied  research  confirmed  the 
basic  weakness  of  co-ordination  design  methods  based  on 
maximizing  bandwidth. 

Grace  and  Potts(7,8)  and  Herman  et  al(’,10),in  and 
around  1964,  gave  a  more  mathematical  treatment  to  Pacey's 
diffusion  theory  and  showed  that  the  model  can  be  described 
by  a  one-dimensional  pseudo-diffusion  equation. 

Hillier  and  RotheryC11),  in  1967,  concluded  that 
platoon  dispersion  can  be  taken  into  account  in  the  design 
of  co-ordination  and  that  the  offset  setting  for  minimum 
delay  at  the  downstream  signal  up  to  305  metres  (1,000  feet) 
is  a  linear  function  of  the  distance  from  the  upstream 
signal . 

Robert son ( 1 2 ) ,  in  1967,  introduced  a  recurrence 
relationship  to  predict  platoon  dispersion  in  his 
contribution  to  the  discussion  on  Area  Control  Of  Road 
Traffic.  This  work  was  carried  out  at  the  Road  Research 
Laboratory  and  is  further  repor ted ( 1 3 , 1  4  ,  1  5  )  in  the 
description  of  the  TRANSYT  (TRAffic  Network  StudY 
Tool)method  of  co-ordinating  traffic  signals.  The  model  uses 
a  recursive  smoothing  equation  and  is  iterative  in  nature. 
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Although  it  appears  that  Robertson  devised  the  method 
empirically  it  has  been  shown  elsewhere ( 1 6 )  that  Robertson’s 
formula  is  similar  to  that  of  Pacey's,  except  that  Pacey's 
transformed,  normal  distribution  of  journey  time  is  replaced 
by  a  geometric  distribution. 

Seddon ( 1 7 , 1 8 , 1 9 )  in  1971/72  attempted  to  validate  the 
three  basic  platoon  prediction  models  proposed  by  Lighthill 
and  Whitham,  Pacey,  and  Robertson,  respectively.  He  found 
that  the  Lighthill' s  and  Whitham' s  theory  could  provide 
reasonable  platoon  prediction  if  an  accurate 
flow-concentration  relationship  is  available.  Pacey's  model 
preformed  well  for  medium  traffic  but  was  not  as  accurate  as 
Lighthill’s  and  Whitham's  in  the  case  of  heavy  traffic 
conditions.  For  best  results,  it  was  concluded  that  all 
models  should  be  calibrated  for  field  conditions. 

Hartley  and  Powner ( 2  0 ) ,  and  Tracz(21),  in  1971, 
suggested  that  a  rectangular  distribution  of  journey  time 
gives  a  similar  predicted  arrival  pattern  to  the  transformed 
normal  or  geometric  distribution  (i.e.  Pacey  or  Robertson). 
This  is  important  because  it  would  make  the  prediction  of 
platoon  dispersion  even  faster  than  by  use  of  the  Pacey  or 
Robertson  model. 

Lam(22),  in  1977,  presented  the  results  of  experimental 
studies  on  the  accuracy  and  sensitivity  of  Robertson's 
model,  and  also  its  calibration  against  actual  field  data  in 
Toronto.  The  results  indicated  that  platoon  data  predicted 
by  Robertson's  model  provided  a  reasonable  fit  with  the 
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observed  data,  but  the  accuracy  could  be  improved  by  the 
introduction  of  a  special  equation  to  initialize  the  platoon 
prediction  process.  Without:  the  initialization  equation, 
there  could  be  significant  errors  in  the  predicated  platoon 
leading  edges,  thus  leading  to  sub-optimal  co-ordination 
designs . 

Among  the  various  studies  in  platoon  dispersion,  the 
research  by  Lighthill  and  Whi tham ( 1 , 2 ) ,  Pacey(3)  and 
Robert  son ( 1 2 )  represent  major  advances  in  the  development  of 
platoon  dispersion  models.  The  work  by  Robertson  is  by  far 
the  most  significant  in  terms  of  practical  application.  This 
is  indicated  by  the  incorporation  of  his  platoon  dispersion 
model  in  such  signal  optimization  programs  as  TRANSYT  (,3) 
and  the  GLC  Combination  Method(23).  The  effectiveness  of 
these  programs  for  both  conventional  and  computerized  signal 
control  systems  has  been  demonstrated  in  a  number  of  field 
appl icat ions ( 2 4 , 2 5 ) ,  and  this  effectiveness  is  partially  due 
to  the  accuracy  of  Robertson's  platoon  model  in  describing 
traffic  flow  patterns  in  a  signal  network. 

2.2  Theoretical  Aspects  Of  Three  Basic  Platoon  Dispersion 
Models 

2.2.1  Lighthill's  and  Whitham's  Wave  Theory 

Lighthill  and  Whithaml1,2)  put  forth  the  theory  that 
kinematic  waves  exist  in  a  traffic  platoon  (or  hump  as  they 

if  there  exists  a  functional  relationship 


described  it) 


16 


between  the  following  explanatory  variables: 


q  =  uk 

where  :  q  =  flow  (vehicles  passing  a  given  point  in 

units  of  time) 

k  =  concentration  (vehicles  present  per  unit 
length  of  road),  and 

u  =  speed  (along  the  road) . 

The  approach  unifies  speed-flow  relationships  at  low 
concentrations  and  headway  phenomena  at  high  concentrations. 
On  the  assumption  of  a  functional  relationship,  the  wave 
theory  follows  from  the  equation  of  continuity  alone. 
Accordingly,  the  waves  are  described  as  "kinematic"  (pure 
motion  without  reference  to  the  masses  or  forces  involved) 
as  opposed  to  "dynamic"  (motion  and  the  equilibrium  of 
systems  under  the  action  of  forces)  waves  which  depend  on 
Newton's  second  law  of  motion  (mass  times  acceration  equals 
force )  . 

Several  characteristics  of  the  flow-concentration  (q  - 
k)  curve,  as  defined  by  Lighthill  and  Whitham,  are 
illustrated  graphically  in  Figure  2.1  and  are  summarized 
below : 

1.  For  any  point  on  the  flow-concentration  curve,  the 
radius  vector  represents  the  traffic  speed  Cus ) . 


2.  The  tangent  to  this  point  on  the  flow-concentration 
curve  represents  the  wave  velocity  (uw) .  It  is  given  by 
differential  of  q  with  respect  to  k  (dq/dk)  at  that 
point.  This  slope  is  smaller  than  the  traffic  speed 
(TFS  )  ,  provided  that  the  traffic  speed  decreases  with  an 
increase  in  concentration.  This  is  acceptable  for 
crowded  roads  and  means  that  disturbances  are  propagated 
backwards  to  the  stream  of  vehicles. 


I 


(  SCALE  OF  ORDER  30  KILO  METERS  ) 
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CONCENTRATION  TIME 

ANALYSIS  OF  TRAFFIC  SHOCKWAVE  USING  A  FLOW- CONCENTRATION  CURVE  AND 
A  TIME- DISTANCE  DIAGRAM.  THE  WAVE  SPEEDS  IN  THE  Tl ME  -  DISTANCE  DIAGRAM 
ARE  DRAWN  PARALLEL  TO  THE  RESPECTIVE  TANGENTS  (JULw  0  Pt.  B  AND  JLLwQ 
Pt.  A  )  OF  THE  FLOW- CONCENTRATION  DIAGRAM. 


PROGRESS  OF  A  TRAFFIC  HUMP  WITH  TIME 


FIGURE  2 . 1 

LIGHTHILL’S  AND  WHITHAM'S  WAVE  THEORY 
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Kinematic  waves  do  not  disperse  as  other  waves  do,  but 
they  suffer  a  change  in  form  due  to  the  dependence  of  the 
wave  speed  on  the  flow  carried  by  the  waves.  Accordingly, 
wave  forms  may  develop  discontinuities  due  to  the  overtaking 
of  slower  waves  by  faster  ones  (with  lower  flow  values). 
These  were  described  by  Lighthill  and  Whitham  as  "shock 
waves"  since  they  were  formed  in  the  same  way  as  that  of 
shock  waves  in  gases.  For  example,  on  the  q  -  k  curve  in 
Figure  2.1,  point  A  represents  a  condition  where  traffic  is 
flowing  near  capacity  and  the  speed  is  reduced  to  a  value 
well  below  the  free-flow  speed.  Point  B,  on  the  other  hand, 
represents  a  condition  where  traffic  flows  at  a  somewhat 
higher  speed  because  of  the  lower  density.  Tangents  at 
points  A  and  B  represents  the  wave  velocities  for  these  two 
conditions.  Now,  if  the  faster  flow  of  point  B  occurs  later 
in  time  than  that  of  point  A,  the  waves  of  point  B  will 
eventually  catch  up  with  those  of  point  A.  This  is 
illustrated  in  the  time-distance  diagram  of  Figure  2.1.  The 
intersection  of  these  two  sets  of  waves  has  a  slope  equal  to 
the  chord  connecting  the  two  points  on  the  q  -  k  curve,  and 
this  intersection  represents  the  path  of  the  shock  wave. 
Note  that  the  velocity  of  the  shock  wave  is  often  negative 
with  respect  to  the  roadway  and  is  always  negative  with 
respect  to  the  traffic. 

The  law  of  motion  of  kinematic  shock  waves  is  derived 
from  considerations  of  conservation  of  vehicles,  as  was  the 
law  of  continuous  waves.  If  the  flow  and  concentration 
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values  (qi,  k^  on  one  side  of  the  shock  waves,  and  the  flow 
and  concentration  (q2,  k2)  on  the  other  side  of  the  shock 
wave,  both  which  move  with  speed  (v),  then  the  quantity 
crossing  per  unit  time  may  be  written  as  : 

q i  "  vk,  =  q 2  -  vk2 

This  gives  the  velocity  of  the  shock  wave  as  : 

v  =  (q2  -  q i )  /  ( k 2  -  k, ) 

The  velocity  of  the  shock  wave,  by  definition,  is  the  slope 
of  the  chord  joining  the  two  points  on  the  q  -  k  curve  (for 
a  given  value  of  x),  which  corresponds  to  the  states  ahead 
of,  and  behind  the  shock  wave  when  it  reaches  x. 

The  waves  on  the  time-distance  diagram  as  presented  are 
not  trajectories  of  vehicles,  but  lines  of  constant  flow  and 
thus  lines  of  constant  speed.  The  individual  vehicles  have  a 
greater  velocity  than  the  waves,  because  the  speed  of 
vehicles  stream  is  represented  by  the  radius  vector,  whereas 
the  velocity  of  the  waves  is  represented  by  the  tangent. 

Lighthill  and  Whitham  applied  such  analysis  techniques 
to  the  progress  of  a  traffic  "hump",  a  hump  being  a  parcel 
of  increased  density,  such  as  might  occur  on  a  freeway, 
flowing  at  a  constant  level  along  the  main  stream,  where 
there  is  a  short-term  influx  of  substantial  proportions  at 
an  on-ramp.  Applying  such  analysis,  traffic  waves  associated 
with  the  formation  of  a  hump  could  be  determined.  The  speed 
of  the  front  of  the  hump  can  be  obtained  immediately  from 
Figure  2.1,  as  the  velocity  of  the  earliest  wave 
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representing  the  increased  flow.  The  construction  of  the 
path  of  the  shock  wave,  as  discussed  earlier,  is  also  shown 
in  Figure  2.1. 

In  summary,  to  obtain  the  kinematic  wave  and  shock  wave 
speeds,  a  detailed  knowledge  of  the  q  -  k  relationship  for 
the  roadway  in  question  is  necessary,  and  this  relationship 
has  to  be  assumed  constant  over  a  given  length  of  road  at  a 
given  time.  This  assumption  would  seem  reasonable  provided 
no  change  in  the  composition  of  the  traffic  stream  (i.e.  the 
proportion  of  slower  moving  trucks  and  buses)  takes  place 
during  the  period  considered.  Clearly  this  composition  and 
the  motivation  of  drivers  is  likely  to  change  during  the 
day,  and  hence  different  q  -  k  relations  would  be  needed. 
However,  during  a  given  peak  hour,  it  seems  reasonable  to 
assume  that  it  would  remain  constant.  From  a  plot  of  wave 
paths  of  both  the  continuous  waves  and  shock  waves, 
flow/time  relationships,  and  thus  the  platoon  structure  at 
any  point  along  a  road  can  be  determined.  It  thus  appears 
that,  provided  the  various  assumptions  are  acceptable,  the 
Lighthill’s  and  Whitham's  wave  theory  offers  a  method  of 
predicting  the  behaviour  of  platoons  as  they  progress  along 
the  road  from  traffic  signals. 

The  following  are  the  three  major  criticisms  of  the 
Lighthill's  and  Whitham's  kinematic  wave  theory: 


1.  The  theory  requires  an  accurate  model  of  q  -  k 
relationships  which  is  difficult  to  obtain,  especally 
for  urban  arterials  where  stop-start  conditions  prevail. 
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2.  The  predicted  platoon  does  not  increase  in  length  as 
such  a  platoon  does  in  reality,  due  to  a  weakness  in  the 
theory . 

3.  The  plotting  of  the  various  wave  paths  is  a  very  time 
consuming  process. 

As  far  as  the  Researcher  knows,  the  Lighthill's  and 
Whitham’s  wave  theory  has  not  been  used  in  any  signal 
control  design  procedure  to  date. 

2.2.2  Pacey's  Platoon  Diffusion  Theory 

Pacey(3)  assumed  that  there  is  no  interference  with 
overtaking  (passing  is  possible  at  will),  and  that  any  given 
vehicle  proceeds  with  the  same  velocity  irrespective  of  the 
number  or  distribution  of  vehicles  on  the  road.  This  is 
likely  to  be  the  case  where  flows  are  light  on  wide  roads, 
but  it  is  unlikely  that  free  overtaking  will  be  possible  in 
congested  urban  conditions. 

Assuming  the  distribution  of  velocities  of  vehicles  to 
be  normal,  it  is  possible  to  evaluate  the  distribution  of 
journey  times  between  the  two  observation  points.  Pacey’s 
platoon  diffusion  may  be  descibed  as  follows:  Let 

f(u)du  =  the  probability  distribution  of  platoon 

speeds,  assumed  to  be  a  normal  distribution, 

r  =  travel  time  between  the  two  observation 
points , 

=  D/u  where  D  =  distance  between  the  two 
observation  points,  and 

g(r)dr  =  distribution  of  travel  times. 
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The  probability  that  a  vehicle  travel  time  lies  between 
r  and  (r  +  dr)  is  the  same  as  the  probability  that  the 
corresponding  speed  lies  between  u  and  (u  -  du) .  Noting  that 
u  =  D  /  r, 

du  =  (  D  /  r 2  )  dr,  and 

g(r)dr  =  f(u)du  =  f[(D/r)  (D/r2)  dr]. 

Now,  because  speeds  are  assumed  to  be  normally  distributed, 

f(u)  =  1/(a  i/(27r))  exp[  -  (  u-u)  2/2cj  2  ]  , 

where  o  is  the  population  standard  deviation  of  speed. 
Therefore , 

g(r)  =  D/a  r2  /(27r)  exp[-(D/r  -  D/r)2/2a2]. 

Letting  s  =  a  /  D, 

g(r)  =  1/r2  s  V{2tt)  exp{[-(1/r  -  u/D)2]/2s2}. 

Cons idera t i ng  now  the  two  observation  points,  the 
number  of  cars  passing  the  first  point  is  q,dt  in  the 
interval  [  t,  (t  +  dt)  ].  Of  this  flow,  the  number  of  cars 
passing  the  second  observation  point  will  be  q , ( t ) g ( r ) dtdT 
in  the  interval  [  (t  +  r),  (t  +  r  +  dT)  ]. 

The  total  number  of  cars  passing  point  two  in  the 
interval  [T,  (T  +  dT)]  is  q2(T)dT  =  / tq , ( t )g (T-t )dtdT,  the 
integral  being  overall  values  of  t  for  which  qi(t)  exceeds 
zero.  As  the  flow  pattern  will  not  be  represented  by  a 
continous  curve  whose  equation  is  known  but  by  a  histogram, 
it  is  more  convenient  to  use  a  corresponding  formula  for  the 
discrete  case.  Thus  the  following  expression  is  employed: 
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q2(  j) 

where  i 
i 

j 

q2  ( j  ) 
q  i  (  i  ) 
g( j-i> 

n 


n 

=  Z  q i ( i )g( j-i ) , 
i  =  1 

=  count  of  the  first  observation  point, 

=  count  of  the  second  observation  point, 

=  flow  past  point  2  during  the  j-th  interval, 

=  flow  past  point  1  during  the  i-th  interval, 

=  probability  of  a  journey  time  of  a 
duration  (j  -  i)intervals,  and 

=  number  of  intervals  considered  (e.g.  number  of 
increments  in  the  cycle  length) . 


Putting  this  formula  into  words,  the  flow  in  the  j-th 
interval  at  the  second  observation  point  equals  the  sum  over 
all  the  values  of  i,  of  the  flow  in  the  i-th  interval  at  the 
first  point  multiplied  by  the  probability  of  a  journey  time 
(j  -  i)  intervals. 

Pacey  found  that  in  general,  the  diffusion  theory  fits 
the  observed  data  fairly  well,  but  that  it  is  slightly 
inferior  to  the  Lightham’s  and  Whitham’s  wave  theory  in 
extreme  traffic  conditions.  For  this  reason,  Pacey's  platoon 
diffusion  theory  has  not  been  used  in  any  signal  control 
design  procedure. 


2.2.3  Robertson's  Recurrence  Model 

Robert son ( 1 2 )  introduced  his  method  of  predicting 
platoon  behaviour,  in  his  contribution  to  the  discussion  at 
the  Institute  of  Civil  Engineers,  Joint  Symposium  On  Area 
Control  Of  Road  Traffic  (London,  England)  in  1967.  His  work 
was  carried  out  at  the  Road  Research  Laboratory  and  is 
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further  reported(13)  in  the  description  of  the  TRANSYT 
(TRAffic  Network  StudY  Tool)  method  of  coordinating  traffic 
signals.  The  method  is  simple  to  apply  and  makes  use  of  an 
empirical  recurrence  relationship. 

This  relationship  can  be  illustrated  by  reference  to 
Figure  2.2.  The  unit  flow  in  the  i-th  interval  of  the 
initial  cyclic  flow  profile!26)  is  multiplied  by  a  factor  F 
to  give  a  portion  of  the  unit  flow  in  (i+/?T)  interval  in  the 
predicted  cyclic  flow  profile.  In  addition,  (i+/?T-1) 
interval  of  the  predicted  cyclic  flow  profile  is  multiplied 
by  a  factor  (1  -  F)  and  added  to  the  (i+/?T)  interval  of  the 
predicted  cycle  flow  profile.  By  repeating  this  process  for 
all  flow  intervals,  the  downstream  arrival  pattern  can  be 
derived.  In  summary,  the  factor  F  controls  the  rate  at  which 
a  platoon  disperses  while  the  effect  of  introducing  j3T  is  to 
offset  the  predicted  pattern  in  time  by  an  equal  amount,  or 
to  ensure  that  no  vehicles  arrive  at  the  downstream  stopline 
until  pT  time  intervals  have  elapsed.  Robertson's  recurrence 
model  takes  the  following  mathematical  form  : 


DIRECTION  OF  TRAFFIC  FLOW 
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ROBERTSON’S  RECURRENCE  MODEL 
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q' 

( i+0T) 

=  Fq  +  ( 1  -  F)q' 

(i)  (i+0T-1) 

where  : 


i  = 

1,2,... n;  where  n  =  number  of  time  intervals 
comprising  the  cycle  length  at  the  upstream 
signal , 

q  = 

i 

the  flow  in  the  i-th  time  interval  of  the 
initial  cyclic  flow  profile  (upstream  flow  in 
interval  i  )  , 

q;  = 

i 

the  flow  in  the  i-th  time  interval  of  the 
predicted  cyclic  flow  profile  (predicted 
downstream  flow  in  interval  i), 

T  = 

the  average  journey  time  over  the  distance 
for  which  the  platoon  dispersion  is  being 
calculated  (meaured  in  the  time  interval 
used  for  q  ) , 
i 

P  = 

(beta),  an  empirical  travel  time  factor 
expressed  as  the  ratio  between  the  platoon 
leader  travel  time  and  the  average  travel 
time  of  the  entire  platoon  (Robertson 
suggested  a  value  of  0.8), 

F  = 

an  empirical  smoothing  factor  given  by  the 
following  expression, 

F  = 

1 

1  +  a/3T 

a  = 

(alpha),  an  empirical  dispersion 
factor  accounting  for  the  degree 
of  platoon  dispersion  between 
signals  (Robertson  suggested  an 
value  of  0.50), 

K  = 

=  a  *  P  (Robertson's  Standard 

-  K  value  =  40) 

100 

j3T  =  integer  value. 
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For  no  dispersion  of  traffic  platoons,  alpha  would 
assume  a  value  of  zero,  beta  would  assume  a  value  of  one, 
and  the  smoothing  factor  F  would  equate  to  a  value  of  one. 

In  the  limit,  as  the  number  of  intervals  approaches 
infinity,  the  predicted  cycle  flow  profile  approaches  zero, 
but  since  the  time  period  is  not  infinity  but  rather  a  cycle 
length,  small  errors  in  the  predictive  process  are 
introduced.  It  should  also  be  noted  that  Robertson’s 
equation  requires  the  flow  value  in  the  previous  downstream 
interval  to  predict  the  value  in  the  present  interval.  This 
information  is  not  available  in  computing  the  first 
downstream  interval.  Hence  Robertson's  equation  will 
consistently  underestimate  the  downstream  flow  rates. 

It  is  important  to  note  that  the  exponentially 
smoothing  process  described  above  relates  to  the  average 
platoon  behaviour.  In  reality,  individual  platoons  will  vary 
randomly  around  the  mean.  (It  should  be  pointed  out  that  the 
equation  on  page  141  of  Reference  (12)  has  been  misprinted.) 

Robert  son ( 1 3 )  stated  that  it  would  be  reasonable  to 
expect  that  the  value  of  F  would  be  a  function  of  various 
site  factors  such  as  road  width,  gradient,  parking  acitiviy, 
opposing  flow  level,  traffic  composition,  etc.,  rather  than 
simply  the  average  journey  time. 

Although  it  appears  that  Robertson  devised  the  method 
emperically,  its  undoubted  success  in  signal  control  design 
procedures  suggests  that  it  must  have  a  sound  theoretical 
basis.  Seddon(19)  has  shown  that  Robertson's  formula  is  the 


. 
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same  as  that  of  Pacey,  except  that  Pacey's  transformed 
normal  distribution  of  journey  g(j-i)  is  replaced  by  a 
probability  function  F(1  -  F)(j_t‘i;.  This  latter  function 
is  a  geometric  distribution  which  is  commonly  used  for 
representing  the  number  of  failures  in  a  series  of  trials  to 
the  first  success.  In  this  case,  it  is  the  probability  that 
a  vehicle  passing  the  first  point  in  the  i-th  interval  will 
pass  the  second  point  in  the  j-th  interval.  Seddon  therefore 
concludes  that  the  success  of  Robertson's  Recurrence  Model 
can  be  attributed  to  the  approximation  of  the  geometric 
distribution  to  the  true  distribution  of  journey  times. 
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2.3  Mathematical  Extension  To  Robertson’s  Recurrence  Model 

As  outlined  in  Section  2.2.3,  Robertson's  Recurrence 
Model  for  platoon  dispersion  makes  use  of  the  following 
simple  recurrence  relationship: 

q’  =  Fq  +  ( 1 -F ) q ’  ,  (  1  ) 

(i+0T)  (i)  (i+0t-1) 

where  : 


i  = 

1,2,...n;  where  n=number  of  time  intervals 
comprising  the  cycle  length  at  the  upstream 
signal , 

q 

i 

the  flow  in  the  i-th  time  interval  of  the 
initial  cyclic  flow  profile  (upstream  flow 
in  interval  i ) , 

q’  = 
i 

the  flow  in  the  i-th  time  interval  of  the 
predicted  cyclic  flow  profile  (predicted 
downstream  flow  in  interval  i), 

T 

the  average  journey  time  over  the  distance 
for  which  the  platoon  dispersion  is  being 
calculated  (measured  in  the  time  intervals 
used  for  q  ) , 
i 

& 

(beta),  an  empirical  travel  time  factor 
expressed  as  the  ratio  between  the  platoon 
leader  travel  time  and  the  average  travel 
time  of  the  entire  platoon  (Robertson  sug¬ 
gested  a  value  of  0.8), 

F 

1 

1  +  aj3T 

a  = 

(alpha),  an  empirical  dispersion  factor 
accounting  for  the  degree  of  platoon  dis¬ 
persion  between  signals  (Robertson  sug¬ 
gested  an  value  of  0.50), 

K 

a*(3  (Robertson’s  Standard  K  value  =  40),  and 

100 

0T  = 

integer  value. 
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Robertson's  equation  requires  the  flow  value  in  the 

previous  downstream  interval  [  q *  ]  to  predict  the 

(  i  +  /3T-  1  ) 

value  in  the  oresent  interval  [  q ’  ].  This  information 

( i  +  /?T) 

is  not  available  in  computing  the  flow  during  the  first 
iteration  and  is  assumed  to  take  a  zero  value.  Therefore 
Equation(l)  takes  the  following  form  for  the  first  interval 
only : 

q'  =  Fq  (2) 

(1+0T)  (1) 

and 

qT  =  Fq  +  ( 1-F)qT  (3) 

(i+£T)  (i)  (i+£T-1) 

for  i=2 , 3 , 4 , . . . ,n  (for  all  succeeding  intervals). 

The  use  of  equation  (2)  underestimates  the  flow  rate  during 
the  first  interval,  and  this  effect  is  perpetuated 
throughout  the  entire  prediction  process  due  to  its 
recursive  nature.  To  provide  some  insight  into  this  volume 
error,  an  expression  is  derived(27)  for  determining  the 
discrepancy  between  initial  and  predicted  cyclic  flow 
prof i les : 

q'  =  Fq  +  ( 1-F)q'  (  1  ) 

(i+0T)  (i)  (i+0T-1) 


By  summing  q'  over  the  entire  cycle  which  consists  of 

(i+j3T) 

n  intervals,  Equation  (1)  can  be  re-written  as  follows  : 


1  I 


- 
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n  n  n 

Z  q'  =  F  Z  q  +  ( 1-F)  Z  q?  (4) 

i  =  1  (i+0T)  i  =  1  (i)  i=  1  (i+0T-1) 


n 

The  expansion  of  the  term  Z  qr  is  as  follows  : 

i  =  1  (  i  +  /3T-  1  ) 


n 

Z  q'  =  q'  +  q’  +  q'  + 

i=  1  (  i  +  /3T-  1  )  (pT)  (1+PT)  (2  +  fiT) 


q'  +  . . .  +  q'  .  (5) 

(3+£T)  (n+jBT-1) 


Adding  q'  -  qf  which  equals  to  a  zero  value  to 

(n+pT)  (n+PT) 

Equation  (5)  we  obtain  the  following  : 


n 

Z  q '  =  q '  +  q  ’  +  q  ’  + 

i=  1  (i  +  PT-1)  (P  T)  (1+PT)  (2+fiT) 


qT 

(2+J3T) 


+  q' 

(n+/3T~1  ) 


q'  "  q' 

(n+PT )  (n+PT) 


(6) 


n 

=  q'  +  Z  q'  -  q’  .  (7) 

(P T)  i= 1  (i+0T)  (n+pT) 


£  .  . 
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By  replacing  the  equality  defined  in  Equation  (7)  into 
Equation  (4)  we  obtain  the  following  : 


n  n 

Z  qf  =  F  Z  q  +( 1-F) [  q'  + 

i= 1  (i+0T)  i = 1  (i)  (0T) 


n 

Z  q'  -  q'  ] 

i  =  1  (i+0T)  ( n  +  /3T ) 


n 

=  F  Z  q  +( 1-F) [  q'  -  q?  ]  + 

i= 1  (i)  (0T)  (n+0T) 


n 

( 1-F)  Z  q' 

i  =  1  (i+0T) 


n  n 

Z  q'  -(1-F)  Z  q’ 

i  =  1  (i+0T)  i  =  1  (i  +  j3T) 


n 

=  F  Z  q  +( 1-F) [  q'  -  q'  ] 

i=  1  (i)  (^T)  (n  +  /3T) 


n  n 

F  Z  q'  =  F  Z  q  +(1-F)[  qT  -  q’  ]  (8) 

i  =  1  ( i  +  /3T )  i  =  1  (i)  (0T)  (n  +  /3T) 


Dividing  by  the  factor  F  in  Equation  (8) 
following  : 


we  obtain  the 
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n  n 

I  q'  =  I  q  +((1-F)/F)[  q'  -  qf  ]  (9) 

i=  1  ( i  +/3T )  i=  1  (i)  (j3T)  (n+0T) 


Under  steady  state  conditions,  if  it  exists,  the  value 

of  q'  and  q'  are  identical  and  therefore  the  last 

'(0 T)  _  (n+PT) 

term  in  equation(9)  is  reduced  to  zero.  Under  this  traffic 

condition  there  will  be  no  volume  error  between  the  initial 

n 

cyclic  flow  profile  [  Z  q  ]  and  the  predicted  cyclic 

n  i= 1  (  i  ) 

flow  profile  [  Z  q'  ].  In  practice,  however, 

i = 1  (i+PT) 

Robertson's  equation  assumes  that  q’  has  an  assumed 

( pT ) 

value  of  zero,  with  the  value  of  q'  being  calculated 

( n+pT) 

in  the  final  step  of  the  iteration  process.  This  difference 

in  value  therefore  leads  to  a  finite  volume  error. 

To  understand  the  potential  impact  of  this  finite 

volume  error  let's  consider  the  feasible  values  that  (1-F)/F 

and  [  q*  -  q'  ]  may  assume  under  realistic  traffic 

(P T)  (n+0T) 

flow  conditions.  For  a  typical  lane  lets  assume  that  the 
smoothing  factor  F  takes  on  a  value  of  0.20  and  q' 

( n  +  j3T ) 

takes  on  a  value  of  1.70.  The  value  of  1.70  is  based  on  a  2 
second  time  slice,  saturation  flow  rate  =  1800  vehicles  per 
hour  green  and  85%  degree  of  saturation.  (Saturation  flow  is 
defined  as  the  number  of  passengers  car  units  which  can 
discharge  at  a  stop-line  of  an  "ideal"  intersection  approach 
(width  3.0  to  3.5m),  moving  straight  ahead  ( i .  e .  ,  no  turning 
movements  are  involved),  with  no  additional  traffic  friction 
(i.e.,  no  bus  stops)  under  the  given  set  of  weather 
conditions,  during  the  given  length  of  the  effective  green 
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interval.  The  concept  of  saturation  flow  assumes  that  after 
the  initial  hesitation  following  the  display  of  green, 
traffic  discharges  at  a  constant  rate,  until  shortly  after 
the  beginning  of  amber,  when  a  sharp  drop  occurs.  Effective 
green  time  is  therefore  defined  as  the  horizontal  dimension 
of  the  rectangular  transformation  of  the  area  encompassed  by 
the  saturation  flow  diagrams.)  Given  this  condition  the 
predicted  downstream  cyclic  flow  profile  is  in  error  by  6.8 
vehicles.  At  a  100  second  cycle  length  (10  seconds 
intergreen  period,  80/20  split  between  major  and  minor 
street)  the  average  volume  of  traffic  per  cycle  would 
approach  to  36.00  vehicles.  However,  at  85%  degree  of 
saturation,  the  predicted  downstream  volume  is  36.00  -  5.40 
=  30.60  vehicles.  The  resulting  degree  of  saturation  based 
on  predicted  volumes  (30.60  -  6.80  =  23.80)  is  calculated  to 
be  66%.  The  potential  impact  of  a  28.6%  error  in  volume 
estimate  due  to  the  prediction  process  can  best  be 
illustrated  by  making  reference  to  the  following  TRANSYT  (V5 
-  V6)  expression  for  random  delay: 

Dr  =  (X) 2 

4 ( 1 -X) 

where:  Dr  =  random  delay  in  units  of  vehicle  hour 

per  hour,  and 

X  =  degree  of  saturation,  defined  as  the 
ratio  of  volume  to  capacity. 
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Hence  for  X=0.66  and  0.85  the  respective  values  of  Dr 
are  0.3  and  1.2.  The  basic  conclusion,  therefore  is  that  a 
28.6%  error  in  volume  estimate  under  high  degrees  of 
saturation  could  increase  the  random  delay  component  by  a 
factor  of  four.  On  the  other  hand,  there  is  a  different  set 
of  other  traffic  conditions  in  which  these  types  of  errors 
due  to  the  prediction  process,  can  be  minimized.  Two 
conditions  have  been  identified.  The  first  condition 
includes  a  larger  smoothing  factor  which  will  result  in  the 
value  of  (1-F)/F  approaching  zero.  The  second  condition 
includes  very  compact  platoon  in  the  middle  of  a  long  cycle 
length  with  no  flows  on  either  side  of  the  platoon.  The 
leading  edge  of  this  particular  downstream  platoon  would 
arrive  no  earlier  than  (i+/3T)  time  units  with  the  cycle 
length  being  of  sufficient  length  that  the  value  of  q' 

(  i  +/3T ) 

approaches  a  zero  value  at  the  end  of  the  calculation. 

Since  traffic  platoons  moving  between  fixed  time 
signals  are  cyclic  in  nature  (the  period  equals  the  cycle 
length),  and  if  steady  state  conditions  exist,  the  flow  for 
an  interval  will  be  the  same  over  all  cycle  lengths.  Based 
on  this  information  an  initialization  equation  is  derived. 

As  was  indicated  previously,  Robertson's  equation 

requires  that  the  flow  value  in  the  previous  downstream 

interval  [  qf  ]  to  predict  the  value  in  the  present 

(i+0T-1) 

interval  [  q'  ].  For  the  first  downstream  interval 

( i+0T) 

[  q'  ]  the  previous  downstream  interval  [  qf  ] 

(1+0T)  (1+0T-1) 

is  unknown  and  is  therefore  assumed  to  be  zero.  Applying  the 
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summation  of  series  techniqiue  and  if  steady  state 

conditions  exist  (27)  Robertson’s  Recurrence  Model  is 

expanded  to  provide  an  initial  estimate  for  the  value  of 

q'  in  the  predicted  cyclic  flow  profile.  The 

(i+0T) 

mathematical  proof  follows. 


Robertson's  model  is  given  by  the  following  form: 

q'  =  F  q  +  (  1  — F )  q '  .  (  1  ) 

(i  +  0T)  i  ( i  +  /3T- 1  ) 

Redefine  equation  (1)  by  looking  at  the  previous  interval. 

q’  =  F  q  +  (  1-F)  q'  (10) 

(i+0T-1)  (i-1)  (i+0T-2) 

Substitute  equation  (10)  into  equation  (1).  Equation  (1) 

takes  the  following  form: 

q'  =F  q  +( 1-F) [F  q  +( 1-F)  q'  ] 

(i+0T)  i  (i-1)  (i+0T-2) 

=F  q  +F(1-F)  q  +(1-F)2  q'  .  (11) 

i  (i-1)  (i  +  /3T-2) 

Redefine  equation  (10)  by  looking  at  the  previous  interval. 

q'  =  F  q'  +  ( 1-F)  q’  (12) 

(i  +  /3T-2)  (i+0T-2)  (i+/3T-3) 


Substitute  equation  (12)  into  equation  (11).  Equation  (11) 
takes  the  following  form: 


q' 

(  i  +  j3T ) 


=F  q  +F (  1  -F )  q 

i  (i-1) 

(1-F)  q’  3 

( i+0T-3) 

F  q  +  F( 1-F)  q 
(i) 


+  (  1-F) 2 [F  q 

(i-2) 


+  F( 1-F) 2  q 


+ 


+ 


(i-1) 


(  i-2 ) 
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( 1-F) 3  q' 


( i+pT- 3 ) 


(13) 


Re-define  equation  (12)  by  looking  at  the  previous  interval. 

q'  =  F  q  +  ( 1-F)  q’  (14) 

(i+PT-3)  (i-3)  (i+0T-4) 


Substitite  equation  (14)  into  equation  (13).  Equation  (13) 
takes  the  following  form: 

q '  =F  q  +F(1-F)  q  +  F(1-F)2  q  + 

(i  +  0T)  (i)  ( i  —  1  )  (i-2) 


( 1-F) 3 [F  q’  +  ( 1-F)  q'  ] 

(i-3)  (i+PT-4) 


=  F  q  +F( 1-F)  q  +F ( 1 -F ) 2  q  + 

i  (  i-  1  )  ( i-2 ) 


F(  1-F) 3  q  +  (  1-F)  4  q' 

(i-3)  ( i+PT-4 ) 


(15) 


It  follows  that  q'  can  be  structured  in  the  following 

(i+0T) 

general  form: 

q’  =F  q  +F ( 1 -F ) 1  q  +F(1-F)2  q 

(i  +  ^T)  i  ( i - 1 )  (i-2) 


+...+  F ( 1 -F ) n “ 1  q  +(1-F)n  q ’  (16) 

(  i-n+  1  )  (  i  +  /3T-n  ) 


More  specifically,  equation  (16)  translates  to  the 
following : 


n-  1 

q ’  =  I  F( 1-F) j  q  +( 1-F) n  q'  (17) 

(i  +  /3T)  j  =  0  (i-j)  (i+0T-n) 


■ 
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Equation  (17)  and  equation  (1)  are  identical.  If  the  number 
of  time  slices  within  a  cycle  length  is  n  and  we  assume  that 
the  value  from  the  predicted  cyclic  flow  profile  is  the  same 
at  the  beginning  of  each  cycle  then: 

q '  =  q '  =  q ' 

(i+0T-n)  (i+0T)  (i+0T+n) 

=  q '  etc .  (18) 

( i+£T+2n) 

This  equality  follows  from  the  key  assumuption  that 
steady  state  conditions  prevail.  For  interval  i=1  equation 
(17)  and  equation  (18)  become  the  following,  respectively: 


q’ 

( 1+0T) 


n-  1 

Z  F(  1-F) j  q  +( 1-F) n  q1 

j  =  0  (  1  -  j  )  (1  +/?T-n  ) 


(19) 


etc . 


q '  =  q '  =  q '  =  q ' 

(1+£T-n)  (1+0T)  (1+£T+n)  (1+£T+2n) 

(20) 


Substitute  q'  in  equation  (19)  with  q'  as 

(1+j8T-n)  (  i  +j3T ) 

defined  by  the  equality  in  equation  (20). 


n-  1 

qf  =  I  F( 1-F) j  q  +  (  1-F) "  q' 

(1+0T)  j  =  0  (1-j)  (  1+/3T) 


(21  ) 


Therefore  : 
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n-  1 

q'  ~  C (  1  -F ) n  q'  ]=  I  F(1-F)j  q  (22) 

(1+/3T)  (  1 +/3T )  j  =  0  (  1  -  j  ) 


q’  ( 1  - (  1  -F ) n ) 

(  1+0T) 


n-  1 

I  F(  1-F) j  q 

j=o  d-j) 


(23) 


n-  1 

q'  =  I  F(1-F)j  q  /[ 1  - ( 1 -F ) n ]  (24) 

(1+/3T)  j  =  0  (1-j) 


Equation  (24)  can  be  used  not  only  for  initial  interval 
prediction  but  also  for  computing  all  intervals  simply  by 
substituting  the  indices  (1+/3T)  and  (1-j)  by  (i+j3T)  and 
(i-j)  respectively.  However,  this  equation  is  more 
cumbersome  to  use  than  Robertson's  equation,  and  after  the 
first  interval  would  not  improve  the  model  accuracy.  The 
incorporation  of  this  equation  into  the  Robertson  Model  will 
in  subsequent  sections  of  this  thesis  be  referred  to  as  the 
Modi f ied  Robertson  Recurrence  Model .  It  is  anticipated  that 
this  improvement  to  the  model  will  significantly  affect  the 
calibration  process.  Since  the  upstream  and  downstream 
cyclic  flow  profiles  will  have  the  same  volume,  the  result 
will  be  better  predictions  with  respect  to  total  delay 
(random  delay  +  uniform  delay),  number  of  vehicle  stops  and 
performance  index. 


3.  DATA  COLLECTION  EFFORTS 


3.1  Data  Collection  Techniques 


The  collection 

of  real-time 

,  traffic  platoon 

data 

was 

achieved  with  the 

aid  of  a 

Micro-computer 

Cyclic 

Flow 

Profile  Monitor!28). 

Travel  time 

information  was 

derived 

by 

collecting  spot  speed  with  the  aid  of  a  radar  speed  gun. 
Technical  details  of  each  of  these  two  tools  are  discussed 
below . 


3.1.1  Platoon  Data 

Research  at  the  University  of  Alberta  is  currently 
underway  to  develop  a  Micro-computer  Cyclic  Flow  Profile 
Monitor.  The  specific  objective  of  this  research,  which  is 
under  the  auspices  of  both  the  Department  of  Electrical 
Engineering  and  the  Department  of  Civil  Engineering,  at  the 
University  of  Alberta,  is  as  follows: 

"  To  Design  And  Build  A  Micro-computer  Based  Local 
Device  To  Monitor  Cyclic  Flow  Profiles  Within  The 
Context  Of  An  Incident  Detection  System.  " 

With  the  support  of  the  City  of  Edmonton, 
Transportation  Department  (formerly  the  Engineering 
Department),  a  basic  system  has  been  designed  and  built  to 
permit  numerous  users  to  collect  several  different  types  of 
traffic  related  information. 

The  system  configuration  of  this  device  makes  use  of  an 
eight  bit  mic ropocessor  chip.  Erasable  Programable  Read  Only 
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Memory  (EPROM)  provides  program  storage,  while  Random  Access 
Memory  (RAM)  is  used  for  data  collection  purposes.  All 
information  collected  in  RAM  is  dumped  onto  a  sixty  (60- 
minute  cassette  tape,  via  a  standard  tape  recorder.  Once  on 
tape,  information  is  then  transmitted  to  the  University  of 
Alberta  Computing  Centre  from  the  City  of  Edmonton  Traffic 
Control  Centre.  Figure  3.1  conceptually  illustrates  a  block 
diagram  of  the  system. 

The  hardware  of  the  system  is  described  below  and  is 
illustrated  conceptually  in  Figure  3.2. 


1.  An  A . C .  Transformer  is  used  to  convert  12  volts  D.C. 
power  from  a  standard  car  battery  to  120  Volts  A.C. 
power  to  energize  the  microprocessor. 


2.  The  entry  of  the  operating  system  commands  for  the 
microprocessor  itself  is  via  a  standard  Keyboard  Panel. 
The  keyboard  is  used  initially  to  start  the  data 
collection  routine,  with  the  basic  input  being  the  cycle 
time  of  the  upstream  signal  and  a  sequence  set  of 
commands  to  set  up  the  work  space.  A  reset  button  on  top 
of  the  microprocessor  terminates  the  data  collection 
effort.  The  keyboard  is  also  used  to  issue  a  set  of 
commands  to  store  data  onto  a  standard  cassette 
recorder . 


3.  The  interface  card,  ribbon  cable,  and  connector  box 
along  with  the  cable  wires  form  the  system  from  which 
the  collection  of  field  data  is  transmitted  to  the 
microprocssor . 


4.  The  length  of  the  cable  wire  essentially  constrain  the 
maximum  possible  distances  between  observations  points. 
Each  survey  location  (shown  as  Location  #1,  #2  and  #3) 
consists  of  two  channels  for  recording  platoon  data. 
Typically  Location  #1  uses  channels  0  and  1,  Location  #2 
uses  channels  2  and  3  and  Location  #3  uses  channels  4 
and  5.  For  consistency,  the  lower  channel  at  each 
location  usually  represents  the  curb  lane  of  the  roadway 
section . 
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FIGURE  3 . 1 

BLOCK  DIAGRAM  OF  THE 

MICRO-COMPUTER  CYCLIC  FLOW  PROFILE  MONITOR 
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HARDWARE  USED  IN  THE  FIELD 
FOR  COLLECTING  REAL-TIME  TRAFFIC  PLATOON  DATA 


' 
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5.  The  collection  of  the  data  requires  a  surveyor  at  each 
location.  Basic  instructions  are  to  press  the  respective 
channel  number  button(s)  (event  recorder)  each  time  a 
vehicle  passes  the  survey  location. 

After  the  survey,  the  collected  data  was  stored  on  a 
cassette  tape,  and  the  information  was  then  transmitted  to 
the  University  of  Alberta  Computer  from  Edmonton’s  Traffic 
Control  Centre.  To  translate  this  raw  data  into  a  meaning 
full  format  a  computing  program  ( BRAD#+*PLOTLI B )  was  written 
to  summarize  the  data  in  a  tabular  and  graphical  manner.  The 
tabular  summarizes  are  in  the  form  of  an  arrival  matrix 
while  the  graphical  plots  are  in  the  form  of  average  and 
individual  cyclic  flow  profiles. 

Table  3.1  summarizes  a  part  of  a  typical  tabular 
summary  of  measured  traffic  platoon  data.  Characteristics  of 
the  arrival  matrix  are  described  below: 

1.  The  arrival  table  is  prefaced  with  four  header  cards  to 
organize  the  data  base  with  respect  to  location,  survey 
date/time,  run  identifier  and  channel  number. 

2.  Rows  represent  cycles  of  data  collected. 

3.  Columns  represent  successive  time  intervals  (user 
specified)  up  to  the  cycle  length. 

4.  Each  entry  is  the  number  of  vehicles  arriving  in  a  given 
time  interval. 


5.  The  observed  total  number  of  vehicles  recorded  in  each 
time  slice  is  summarized  along  with  an  unbiased  estimate 
of  the  standard  deviation  (The  information  is  the  basic 
input  into  two  computer  programs  which  will  be  discussed 
in  Chapter  4 ) . 


6.  The  totals  for  each  time  interval  are  also  converted 


LOCATION  104  AVENUE  -  116  STREET  ;  EASTBOUNO  VEHICULAR  FLOWS 

SURVEY  DATE/TIME  :  THURSDAY  ,  FEBRUARY  4  ,  1982  ;  7:30  A.M 

RUN  IDENTIFIER  WINTER  CONDITIONS  .  EXHAUST  FUMES  IN  THE  AIR  ,  ICY  PAVEMENT 

CHANNEL  NUMBER  :0  BOTH  LANES  (CHANNEL  0+1  )  80  FEET  DOWNSTREAM  OF  STOPBAR 
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TABLE  3 . 1 

TYPICAL  ARRIVAL  MATRIX 


■ 
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into  an  average  one  lane  flow  rate  (vehicles  per  hour 
green)  along  with  an  unbiased  estimate  of  the  standard 
deviat ion . 

The  graphical  output  from  this  program  includes  the 
following : 


1.  A  plot  of  the  average  lane  flow  rates  with  an  asterisk 
(*)  marker  to  illustrate  the  degree  of  variability  (  + 

or  -  one  standard  deviation)  in  the  data  base.  A  typical 
graphical  summary  of  the  arrival  matrix  is  shown  in 
Figure  3.3 


2.  A  plot(s)  for  a  specified  cycle(s)  to  observe  the  actual 
lane  flow  rates  for  that  cycle  in  comparision  to  the 
average  lane  flow  rates  for  all  cycles  collected.  The 
average  lane  flow  rates  for  all  cycles  collected  is 
indicated  with  an  asterisk  (*)  marker.  A  detail 
graphical  summary  of  the  arrival  matrix  is  shown  in 
Figure  3.4 


FLOW  RATE  -  VEHICLES  PER  HOUR  GREEN 
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FLOW  RfiTE  VERSUS  CYCLE  TIME  (  AVERAGE  CYCLE  ) 

(  *  INDICATES  RANGE  OF  1  STANDARD  DEVIATION  ) 

LOCATION  i  104  AVENUE  -  116  STREET  J  EASTBOUND  VEHICULAR  FLOWS 

SURVEY  DATE/TIME  :  THURSDAY  .  FEBRUARY  4  .  1982  I  7«30  A.M. 

RUN  IDENTIFIER  t  WINTER  CONDITIONS  .  EXHAUST  FUMES  IN  THE  AIR  .  ICY  PAVEMENT 

CHANNEL  NUMBER  fO  BOTH  LANES  (CHANNEL  0*1  )  80  FEET  DOWNSTREAM  OF  STOPBAR 


CYCLE  TIME  -  SECONDS  FROM  START  OF  MAIN  STREET  GREEN 


FIGURE  3.3 


TYPICAL  GRAPHICAL  SUMMARY  OF  THE  ARRIVAL  MATRIX 


FLOW  RATE  -  VEHICLES  PER  HOUR  GREEN  -n  FLOW  RATE  -  VEHICLES  PER  HOUR  GREEN 
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FLOW  RATE  VERSUS  CYCLE  TIME  (  CYCLE  *  2  ) 

l  *  INDICATES  AVERAGE  CYCLE  FOR  CHANNEL  ) 

LOCUTION  i  104  AVENUE  -  116  STREET  »  EASTBOUNO  VEHICULAR  FLOWS 

SURVEY  OATE/TIhE  I  THURSDAY  .  FEBRUARY  4  .  1982  l  7«30  R.H. 

RUN  IDENTIFIER  «  WINTER  CONDITIONS  .  EXHAUST  FUMES  IN  THE  AIR  .  ICY  PAVEMENT 

CHANNEL  NUMBER  lO  BOTH  LANES  (CHANNEL  0»1  )  80  FEET  DOWNSTREAM  OF  ST0P8RR 


CYCLE  TIME  -  SECONDS  FROM  START  OF  MAIN  STREET  GREEN 


W  RATE  VERSUS  CYCLE  TIME  (  CYCLE  «  3  ) 

(  *  INDICATES  AVERAGE  CYCLE  FOR  CHANNEL  ) 

LOCATION  «  104  AVENUE  -  116  STREET  l  EASTBOUNO  VEHICULAR  FLOWS 

SURVEY  DATE/TIME  I  THURSDAY  .  FEBRUARY  4  .  1982  »  7i30  A.M. 

RUN  IDENTIFIER  i  WINTER  CONDITIONS  .  EXHAUST  FUMES  IN  THE  AIR  .  ICY  PAVEMENT 

CHANNEL  NUMBER  lO  BOTH  LANES  (CHANNEL  0*1  )  80  FEET  DOWNSTREAM  OF  STOPBRR 


FIGURE  3.4 


DETAIL  GRAPHICAL  SUMMARY  OF  THE  ARRIVAL  MATRIX 
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3.1.2  Speed/Travel  Time  Data 

To  obtain  an  estimate  of  the  actual  travel  time  between 
two  survey  locations,  the  data  collected  with  spot  speed 
survey  technique  was  translated  to  travel  time  with  the 
uniform  acceleration  model. 

The  standard  radar  gun  which  is  based  on  the 
Doppler-principle  was  used  and  calibrated  prior  to  each  spot 
speed  survey.  A  statistical  analysis  of  each  spot  speed 
survey  date  was  undertaken  using  the  STPK  Package 
(Interactive  Statistical  Graphics  PacKage  (34).  Principal 
outputs  from  this  analysis  included  the  mean  speed 
(kilometres  per  hour),  variance,  and  a  fitted  normal 
probability  density  funtion.  For  all  spot  speed  data 
collected  and  analyzed,  a  Normal  Distribution  was 
statistically  significant  at  the  95%  confidence  level. 
Figure  3.5  illustrates  a  typical  output  from  the  STPK 
analy s i s . 

The  uniform  acceleration  model  is  based  on  the 
assumption  that  acceleration  is  constant.  This  results  in 
the  following  fundamental  relationships  : 


* 
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FIGURE  3.5 

STPK  STATISTICAL  ANALYSIS  OF  TYPICAL  SPOT  SPEED  SURVEY  DATA 
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dv 

—  =  a 

dt 

V 

I  dv  =  jo  a  dt 

Vo 

v  =  v0  +  at 

Jo  dx  =  Jo  ( v o  +  at)  dt 
x  =  v0t  +  1/2  at2 

where  a  =  acceleration 
v  =  speed 

v0  =  initial  speed 
x  =  distance 
t  =  time 

Equation  (1)  to  (3)  are  the  acceleration-time,  speed-time, 
and  di stance- t ime  relationships  respectively  for  uniformly 
accelerated  motion.  The  expression  for  distance  as  a 
function  of  speed  may  be  obtained  from  (1)  as  follows  : 
dv  dx  dv 

dx  dt  dx 

/ o  a  dx  =  /  v  dv 
1 

x  =  — (v2  -  v o )  (4) 

2a 

These  relationships  are  illustrated  in  Figure  3.6 

Based  on  several  spot  speed  survey  mean  speeds  along 
the  roadway  section,  and  the  uniform  acceleration  model,  an 
estimate  of  the  travel  time  between  two  survey  locations  was 


(  1  ) 

(2) 

(3) 


determined . 


Accelerotion  dv/dt  Speed  v  Distonce 
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Speed  v 


dv/dt  -  a 


Time  t 


1 

1 

1 

1 
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1 

dv/dt  -  a 

1 

1 

1 

1 

1 

1 

1 

1 

A 

Speed  v 


FIGURE  3.6 


UNIFORM  ACCELERATION  MODEL 
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3.2  The  Survey  Sites  And  Their  Characteristics 

Table  3.2  provides  an  overview  of  the  survey  sites 
during  the  A.M.  peak  hour  from  the  perspective  of  weather 
conditions,  dispersion  distance,  number  of  platoons 
collected,  total  volume  observed,  and  the  degree  of 
saturation  (ratio  of  volume  to  capacity)  observed.  The 
definition  of  the  three  states  of  weather  conditions  is 
consistent  with  the  research  done  in  Alberta  on  saturation 
flow  ( 3  5 ) . 

"Summer"  driving  conditions  were  defined  as  : 
dry  weather  (i.e.  no  rain  or  snow) 
temperature  form  about  -10°C  to  well  above  0°C 
dry  pavement  (no  snow  residues) 
good  visibility 

"Winter"  driving  conditions  were  characterized  as  : 
dry  weather  (no  snowfall) 
clear  or  cloudy  sky 

temperatures  from  about  -10°C  to  -30°C 

dry  pavement  with  some  packed  snow  residues  but  well 
sanded 

no  snowbanks  at  intersections 
good  visiblity 

"Severe  winter"  conditions  were  characterized  as  : 
clear  or  cloudy  sky 

very  low  temperature  (lower  than  ~30°C) 
visibility  obstructed  by  exhaust  fumes 
or  by : 


, 


t 


TABLE  3.2 


OVERVIEW  OF  THE  SURVEY  SITES  DURING  THE  A.M.  PEAK  HOUR 


Locat ion 

Weather 
Condi t ions 

D i sper s i on 

D i stance 
( meters ) 

Number  Of 

P 1 atoons 
Col lected 

Total 

Vo  1 ume 
( veh i c 1 es 

Degree  Of 
Saturat ion 

)  ( % ) 

104  Av.  - 
116  St . 

Severe 

Winter 

i  06 

21 

860 

92 

104  Av.  - 
116  St . 

Wi nter 

94 

28 

1248 

88 

104  Av.  - 
116  St  . 

( Week 

No.  1  ) 

Summer 

106 

27 

1318 

96 

104  Av.  - 
116  St . 

( Week 

No.  2) 

Summer 

106 

27 

1260 

92 

111  Av.  - 
116  St . 

Summer 

106 

25 

386 

44 

102  A  Av. 
97  St. 

-  Summer 

84 

23 

157 

29 
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heavy  snowfall  or  shortly  after,  before  snow  has  been 
cleared 

temperatures  below  0°C. 

As  can  be  seen,  the  definitions  do  not  provide  for  a 
clear  seasonal  distinction.  "Summer"  driving  conditions  may 
occur  during  the  winter  season  while  "winter"  or  "severe 
winter"  conditions  may  be  experienced  in  the  spring  or  fall. 
Generally,  however,  the  definitions  charac ter i zed  typical 
long  lasting  "summer"  and  "winter"  contitions. 

The  dispersion  distance  generally  varies  between  84  - 
106  metres,  a  distance  which  is  generally  representative  of 
the  majority  of  roadway  links  in  the  CBD,  as  illustrated  in 
Figure  3.7.  The  number  of  platoons  collected  reflects  on 
average  30  -  45  minutes  of  field  observations  per  survey 

site.  These  observations  translated  to  a  volume  of  157  to 
1248  vehicles  with  degrees  of  saturation  varying  from  29  - 
92  % .  In  summary,  more  than  150  platoons  of  data  were 
collected  and  5200  individual  vehicles  were  surveyed  during 
these  surveys. 

3.2.1  104  Avenue  -  116  Street  (Severe  Winter  Conditions) 


1.  Organization  of  data  base  .for  future  reference  : 
LOCATION :  104  Avenue  -116  Street;  Eastbound  Vehicular 
flows.  SURVEY  DATA/TIME :  Thursday,  February  4,1981;  7:30 
A.M.,  RUN  IDENTIFIER:  Severe  Winter  Conditions,  Exhaust 
Fumes  In  The  Air,  Icy  Pavement,  CHANNEL  NUMBERS : 
Locations  #1  (Channel  0  +  1  Combined,  24  Metres  T80 

Feet)  Downstream  Of  The  Signal  Stopbar).  Location  #3 
(Channel  4+5  Combined,  130  Metres  (426  Feet) 


» 

' 


ABSOLUTE  NUMBER  OF  ROADWAY  LINKS  FALLING  IN  THE  CLASS  SHOWN 
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FIGURE  3.7 

HISTOGRAM  OF  THE  LENGTH  IN  METRES  BETWEEN  ALL 
SIGNALIZED  INTERSECTIONS  IN  THE  CENTRAL  BUSINESS 
DISTRICT  OF  EDMONTON,  ALBERTA 


■ 
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Downstream  Of  The  Signal  Stopbar). 


2.  Figure  3.8  illustrates  a  scaled  schematic  drawing  of  the 
104  Avenue  -  116  Street  intersection  and  the  roadway 
section  where  these  cyclic  flow  profiles  were  collected. 

3.  Travel  time  between  Location  #1  and  Location  #3  (a 
distance  of  106  Metres)  is  14.04  seconds,  24.64  feet  per 
second  or  16.80  miles  per  hour. 


4.  A  total  of  21  cycles  of  data  were  collected  with  an 
average  volume  per  cycle  of  40.95  vehicles  for  two  lanes 
with  an  unbiased  standard  deviation  of  7.74  vehicles. 


5.  Signal  timing  information  at  the  upstream  signalized 
intersection  (104  Avenue  -  116  Street);  cycle  length=90 
seconds,  main  street  green  (104  Avenue)=50  seconds,  main 
street  intergreen  (3  seconds  amber,  2  seconds  all-red)=5 
seconds,  cross  street  green  (116  street)=30  seconds,  and 
cross  street  intergreen  (3  seconds  amber,  2  seconds 
all-red)=5  seconds.  Interval  size  for  analysis 
purposes=2  seconds=1  step. 


6.  Characteristics  of  the  survey  site  includes  the 

following  :  major  arterial  roadway  system  leading  into 
CBD ,  a  92%  degree  of  saturation  (a_t  all  survey  sites 
this  degree  of  saturation  i s  implied  for  the  upstream 
s i qnal i zed  inter sec t i on ) ,  no  parking/stopping  along  the 
curb,  very  little  side  friction  (no  driveway),  passing 
freedom  to  manoeuvre  is  seriously  limited  due  to  the 
very  high  degree  of  saturation,  bus  volume  is 
approximately  four  buses  per  hour  with  a  dwell  time  in 
the  range  of  10  -  15  seconds,  and  the  movement  of 

traffic  platoons  can  be  characterized  as  unrestricted 
(the  location  of  the  nearest  downstream  signal  was 
approximately  nine  city  blocks). 

7.  Figure  3.9  presents  graphically  the  survey  information 
for  104  Avenue  -  116  Street  (Severe  Winter  Conditions). 


3.2.2  104  Avenue  -  116  Street  (Average  Winter  Conditions, 
Further  Downstream) 


INFORMATION 
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FIGURE  3.8 

A  SCALED  SCHEMATIC  DRAWING  OF  THE  104  AVENUE  -  116  STREET 
INTERSECTION  AND  THE  ROADWAY  SECTION  WHERE  CYCLIC  FLOW 

PROFILES  WERE  COLLECTED 


- 
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FLOW  RATE  VERSUS  CYCLE  TIME  (  AVERAGE  CYCLE  ) 

(  «  INDICATES  RANGE  OF  1  STANDARD  DEVIATION  ) 

LOCATION  i  104  AVENUE  -  116  STREET  l  EASTBOUNO  VEHICULAR  FLOWS 

SURVEY  DATE/TIME  i  THURSDAY  .  FEBRUARY  4  .  1982  l  H30  A.M. 


TRAVEL  TIME  =  14.04 


SECONDS 


FLOW  RATE  VERSUS  CYCLE  TIME  (  AVERAGE  CYCLE  ) 

(  *  INDICATES  RANGE  OF  1  STANDARD  DEVIATION  ) 

LOCATION  i  104  AVENUE  -  116  STREET  »  EASTBOUNO  VEHICULAR  FLOWS 

SURVEY  DATE/TIME  i  THURSDAY  .  FEBRUARY  4  .  1982  I  7i30  A.M. 

RUN  IDENTIFIER  i  WINTER  CONDITIONS  .  EXHAUST  FUMES  IN  THE  AIR  .  ICY  PAVEMENT 

CHANNEL  NUMBER  t4  BOTH  LANES  (CHANNEL  4*6  )  426  FEET  DOWNSTREAM  OF  STOPBRR 


FIGURE  3.9 


GRAPHICAL  SUMMARY  OF  THE  COLLECTED  SURVEY  INFORMATION 
FOR  104  AVENUE  -  116  STREET  (SEVERE  WINTER  CONDITIONS) 
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1.  Organization  of  data  base  for  future  reference  : 
LOCATION:  104  Avenue  -116  Street;  Eastbound  Vehicular 
flows.  SURVEY  DATE/TIME;  Wednesday,  March  10,1981;  7:30 
A.M.,  RUN  IDENTIFIER:  Winter  Month,  Visual  Inspection  Of 
The  Roadway  Surface  Would  Suggest  Typical  Summer 
Conditions  (Dry  Pavement).  CHANNEL  NUMBERS :  Locations  #1 
(Channel  0  +  1  Combined,  130  Metres  (426  Feet) 

Downstream  Of  The  Signal  Stopbar).  Location  #3  (Channel 
4  +5  Combined,  224  Metres  (734  Feet)  Downstream  Of  The 

Signal  Stopbar ) . 


2.  Figure  3.8  illustrates  a  scaled  schematic  drawing  of  the 
104  Avenue  -  116  Street  intersection  and  the  roadway 

section  where  these  cyclic  flow  profiles  were  collected. 


3.  Travel  time  between  Location  #1  and  Location  #3  (a 
distance  of  94  metres)  is  8.24  seconds,  37.40  feet  per 
second  or  25.50  miles  per  hour. 


4.  A  total  of  28  cycles  of  data  were  collected  with  an 
average  volume  per  cycle  of  44.57  vehicles  for  two  lanes 
with  an  unbaised  standard  deviation  of  8.09  vehicles. 


5.  Signal  timing  information  at  the  upstream  signalized 
intersection  (104  Avenue  -  116  Street):  see  Section 

3.3.1. 


6.  Characteristics  of  the  survey  site:  see  Section  3.3.1. 
Degree  of  saturation  during  the  survey  was  88%. 

7.  Figure  3.10  presents  graphically  the  survey  information 
for  104  Avenue  -  116  Street  (Average  Winter  Conditions, 
Further  Downstream) . 


3.2.3  104  Avenue  -  116  Street  (Week  No.  1) 


1.  Organization  of  data  base  for  future  reference  : 
LOCATION:  104  Avenue  -116  Street;  Eastbound  Vehicular 

flows.  SURVEY  DATE/TIME:  Thursday,  April  16,1981;  7:30 

A.M.,  RUN  IDENTIFIER:  Summer  Conditions,  Dry  Pavement, 
CHANNEL  NUMBERS:  Location  #1  (Channel  0  +  1  Combined,  24 
Metres  T80^  Feet)  Downstream  Of  The  Signal  Stopbar). 
Location  #3  (Channel  4+5  Combined,  130  Metres  (426 


* 
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FLOW  RfiTE  VERSUS  CYCLE  TIME  (  AVERAGE  CYCLE  ) 

(  k  INDICATES  RANGE  OF  1  STANDARD  DEVIATION  ) 

LOCUTION  I  104  HVENUE  -  116  STREET  I  EflSTBOUNO  VEHICULAR  FLOWS 

SURVEY  DATE/TIME  i  WEDNESDAY  .  MARCH  10  .  1983  l  7i30  R.M. 

RUN  IDENTIFIER  i  WINTER  MONTH  /  SUMMER  CONDITIONS  .  DRY  PAVEMENT 

CHANNEL  NUMBER  10  BOTH  LANES  (CHANNEL  0*1  I  426  FEET  DOWNSTREAM  OF  STOPBAR 


TRAVEL  TU*E  =  8.24 


SECONDS 


FLOW  RATE  VERSUS  CYCLE  TIME  (  AVERAGE  CYCLE  ) 

l  *  INDICATES  RANGE  OF  1  STANDARD  DEVIATION  ) 

LOCATION  i  104  AVENUE  -  116  STREET  I  ERSTBOUND  VEHICULAR  FL0W6 
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Feet)  Downstream  Of  The  Signal  Stopbar). 


2.  Figure  3.8  illustrates  a  scaled  schematic  drawing  of  the 
104  Avenue  -  116  Street  intersection  and  the  roadway 
section  where  these  cyclic  flow  profiles  were  collected. 

3.  Travel  time  between  Location  #1  and  Location  #3  (a 
distance  of  106  metres)  is  9.87  seconds,  35.05  feet  per 
second  or  23.9  miles  per  hour. 


4.  A  total  of  27  cycles  of  data  were  collected  with  an 
average  volume  per  cycle  of  48.81  vehicles  for  two  lanes 
with  an  unbiased  standard  deviation  of  4.59  vehicles. 


5.  Signal  timing  information  at  the  upstream  signalized 
intersection  (104  Avenue  -  116  Street):  see  Section 

3.3.1. 


6. 


Characteristics  of  the  survey  site:  see  Section  3.3.1. 
Degree  of  Saturation  during  the  survey  was  96%. 


7. 


Figure  3.11  presents  graphically  the  survey  information 
for  104  Avenue  -  116  Street  (Week  No.  1). 


3.2.4  104  Avenue  -  116  Street  (Week  No.  2) 


1.  Organization  of  data  base  for  future  reference  : 

LOCATION:  104  Avenue  -116  Street;  Eastbound  Vehicular 

flows.  SURVEY  DATE/TIME:  Thursday,  April  23,1981;  7:30 

A. M. (Note:  This  survey  was  duplicated  one  week  later), 
RUN  IDENTIFIER:  Summer  Conditions,  Dry  Pavement,  CHANNEL 
NUMBERS:  Location  #1  (Channel  0  +  1  Combined,  24  Metres 

(80  Feet)  Downstream  Of  The  Signal  Stopbar).  Location  #3 
(Channel  4  +5  Combined,  130  Metres  (426  Feet) 

Downstream  Of  The  Signal  Stopbar). 

2.  Figure  3.8  illustrates  a  scaled  schematic  drawing  of  the 

104  Avenue  -  116  Street  intersection  and  the  roadway 

section  where  these  cyclic  flow  profiles  were  collected. 

3.  Travel  time  between  Location  #1  and  Location  #3  (a 
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(  *  INDICATES  RANGE  OF  1  STANDARD  DEVIATION  I 


LOCATION 
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CHANNEL  NUMBER 
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i  SUMMER  CONDITIONS  .  DRY  PAVEMENT 

«0  BOTH  LANES  (CHANNEL  0*1  )  80  FEET  OOWN8TREAM  OF  STOPBAR 


CYCLE  TIME  -  SECONDS  FROM  START  OF  MAIN  STREET  GREEN 


TRAVEL  TIME  = 


9.  87 
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FLOW  RATE  VERSUS  CYCLE  TIME  ;  AVERAGE  CYCLE  ) 

(  «  INDICATES  RANGE  OF  1  STANDARD  DEVIATION  ) 

LOCATION  1  104  AVENUE  -  116  STREET  l  EASTBOUNO  VEHICULAR  FLOWS 

SURVEY  DATE/TIME  i  THURSDAY  .  APRIL  16  ,  1981  I  7i30  A.M. 

RUN  IDENTIFIER  I  SUMMER  CONDITIONS  .  DRY  PAVEMENT 

CHANNEL  NUMBER  «4  BOTH  LANES  (CHANNEL  4*6  )  426  FEET  DOWNSTREAM  OF  STOPBAR 


FIGURE  3.11 


GRAPHICAL  SUMMARY  OF  THE  COLLECTED  SURVEY  INFORMATION 
FOR  104  AVENUE  -  116  STREET  (WEEK  NO.  1) 
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distance  of  106  metres)  is  9.87  seconds,  35.05  feet  per 
second  or  23.9  miles  per  hour. 


4.  A  total  of  27  cycles  of  data  were  collected  with  an 
average  volume  per  cycle  of  46.70  vehicles  for  two  lanes 
with  an  unbiased  standard  deviation  of  7.66  vehicles. 


5.  Signal  timing  information  at  the  upstream  signalized 
intersection  (104  Avenue  -  116  Street):  see  Section 

3.3.1. 


6.  Characteristics  of  the  survey  site:  see  Section  3.3.1. 
Degree  of  Saturation  during  the  survey  was  92%  (decrease 
of  four  percentage  points  from  the  previous  survey). 

7.  Figure  3.12  presents  graphically  the  survey  information 
for  104  Avenue  -  116  Street  (Week  No.  2). 


3.2.5  111  Avenue  -  116  Street 


1 .  Organization  of  data  base  for  future  reference  : 

LOCATION:  111  Avenue  -116  Street;  Eastbound  Vehicular 
flows.  SURVEY  DATE/TIME:  Friday,  May  6,1983;  7:30  A.M., 
RUN  IDENTIFIER:  Summer  Conditions,  Dry  Pavement,  CHANNEL 
NUMBERS:  Location  #1  (Channel  0  +  1  Combined,  24  Metres 

(80  Feet)  Downstream  Of  The  Signal  Stopbar).  Location  #3 
(Channel  4  +5  Combined,  130  Metres  (426  Feet) 

Downstream  Of  The  Signal  Stopbar). 


2.  Figure  3.13  illustrates  a  scaled  schematic  drawing  of 
the  111  Avenue  -  116  Street  intersection  and  the  roadway 
section  where  these  cyclic  flow  profiles  were  collected. 


3.  Travel  time  between  Location  #1  and  Location  #3  (a 

distance  of  106  metres)  is  11.08  seconds,  31.24  feet  per 
second  or  21.30  miles  per  hour. 

4.  A  total  of  25  cycles  of  data  were  collected  with  an 
average  volume  per  cycle  of  15.44  vehicles  for  two  lanes 
with  an  unbiased  standard  deviation  of  4.06  vehicles. 
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A  SCALED  SCHEMATIC  DRAWING  OF  THE  111  AVENUE  -  116  STREET 
INTERSECTION  AND  THE  ROADWAY  SECTION  WHERE  CYCLIC  FLOW 

PROFILES  WHERE  COLLECTED 


- 


■ 


67 


5.  Signal  timing  information  at  the  upstream  signalized 
intersection  (111  Avenue  -  116  Street):  cycle  length=70 
seconds,  main  street  green  (111  Avenue)=37  seconds,  main 
street  intergreen  (3  seconds  amber,  1  second  all-red)=4 
seconds,  cross  street  green  (116  street)=24  seconds,  and 
cross  street  intergreen  (3  seconds  amber,  2  seconds 
all-red)=5  seconds.  Interval  size  for  analysis 
purposes=2  seconds=1  step. 


6.  Charac ter i st ic s  of  the  survey  site  includes  the 

following  :  major  arterial  roadway  system  leading  into 

CBD,  medium  degree  of  traffic  pressure  as  measured  by  a 
44%  degree  of  saturation,  no  parking/stopping  along  the 
curb  (note:  curb  lane  was  not  surveyed),  very  little 
side  friction  (no  driveway),  passing  freedom  is 

generally  not  restricted  due  to  medium  degree  of 
saturation  (significant  amount  of  passing  was  observed), 
bus  volume  approximately  four  buses  per  hour  with  a 
dwell  time  in  the  range  of  10  -  15  seconds,  and 

unrestricted  traffic  platoon  movement  (the  location  of 
the  nearest  downstream  signal  was  approximately  seven 
city  blocks ) . 


7.  Figure  3.14  presents  graphically  the  survey  information 
for  111  Avenue  -  116  Street. 

3.2.6  102  A  Avenue  -  97  Street 


1.  Organization  of  data  base  for  future  reference  : 

LOCATION:  102  A  Avenue  -  97  Street;  Eastbound  Vehicular 
flows.  SURVEY  DATE/TIME:  Friday,  October  14,1983;  7:25 
A.M.,  RUN  IDENTIFIER:  Summer  Conditions,  Dry  Pavement, 
CHANNEL  NUMBERS:  Location  #1  (Channel  0  +  1  Combined,  21 
Metres  (70  Feet)  Downstream  Of  The  Signal  Stopbar). 
Location  #3  (Channel  4  +5  Combined,  105  Metres  (345 

Feet)  Downstream  Of  The  Signal  Stopbar). 

2.  Figure  3.15  illustrates  a  scaled  schematic  drawing  of 

the  102  A  Avenue  -  97  Street  intersection  and  the 

roadway  section  where  these  cyclic  flow  profiles  were 
collected . 


3.  Travel  time  between  Location  #1  and  Location  #3  (a 
distance  of  84  metres)  is  12.76  seconds,  21.56  feet  per 
second  or  14.70  miles  per  hour  (function  of  the 
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A  SCALED  SCHEMATIC  DRAWING  OF  THE  102  A  AVENUE  -  97  STREET 
INTERSECTION  AND  THE  ROADWAY  SECTION  WHERE  CYCLIC  FLOW 

PROFILES  WERE  COLLECTED 
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space-time  diagram  for  southbound  flows  between  102  A 
Avenue  and  102  Avenue). 


A  total  of  23  cycles  of  data  were  collected  with  an 
average  volume  per  cycle  of  6.83  vehicles  for  two  lanes 
with  an  unbiased  standard  deviation  of  3.07  vehicles. 


5.  Signal  timing  information  at  the  upstream  signalized 
intersection  (102  A  Avenue  -  97  Street)  is  presented  in 
Figure  3.16  as  a  space-time  diagram.  In  reviewing  this 
figure  it  is  important  to  note  that,  by  design,  the 
southbound  green  waveband  is  coordinated  for  a  speed 
23.7  kilometres  per  hour  (14.7  miles  per  hour).  Spot 
speed  surveys  undertaken  immediately  downstream  of  102  A 
Avenue,  indicate  that  vehicles  desire  to  achieve  a  speed 
of  approximately  35  kilometres  per  hour.  Because  of  the 
signal  design  the  phenonenon  of  platoon  compression 
appears  to  be  taking  place  (reverse  of  platoon 

dispersion).  Interval  size  for  analysis  purposes=2 
seconds=1  step. 


6.  Characteristics  of  the  survey  site  include  the  following 
:  major  arterial  roadway  system  in  the  CBD,  very  low 
degree  of  traffic  pressure  as  measured  by  a  29%  degree 
of  saturation,  no  parking/stopping  along  the  curb,  very 
little  side  friction,  passing  freedom  to  manoeuvre  is 
not  restricted  at  all  due  to  the  low  degree  of 
saturation,  and  restricted  traffic  platoon  movement  (the 
proximity  of  the  downstream  signal  affects  the 
dispersion  or  alternately  compression  of  traffic 
platoons).  (Non-platoon  to  a  platoon  situation, 
compression  of  flow  into  a  platoon  or  a  re-grouping  of 
vehicles . ) 


7.  Figure  3.17  presents  graphically  the  survey  information 
for  102  A  Avenue  -  97  Street. 
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SPACE-TIME  DIAGRAM  ALONG  97  STREET  (SOUTHBOUND)  BETWEEN 

102  A  AVENUE  AND  102  AVENUE 
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FLOW  RATE  VERSUS  CYCLE  TIME  (  AVERAGE  CYCLE  ) 

(  *  INDICATES  RANGE  OF  1  STANDARD  DEVIATION  ) 
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FIGURE  3.17 
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4.  THE  ANALYTICAL  TOOLS  FOR  CALIBRATION 


To  assist  in  the  manipulation,  handling,  and  management 
of  the  data  collected  with  respect  to  model(s)  calibration 
three  computing  programs  were  written.  They  are  described 
below : 

4.1  The  GORDONV1  Computing  Program 

The  computing  program  GORDONV1  is  based  on  the  program 
ROBERT,  that  was  developed  at  the  Toronto  Traffic  Control 
Centre  (29).  The  version  of  ROBERT  which  was  obtained  was 
found,  upon  vigours  testing  to  have  several  software  errors 
(/3T  did  not  assume  an  integer  value,  prediction  of  the 
downstream  cyclic  flow  profile  was  not  displaced  by  0T  time 
intervals  and  the  delay/offset  algorithm  was  not 
functioning).  Based  on  the  skeleton  program  of  the  ROBERT 
program,  revisions  where  made  to  suit  the  specific  needs  of 
this  research  and  the  program  was  renamed  G0RD0NV1 . 

The  GORDONV1  computing  program  calculates  the 
downstream  cyclic  flow  profile  from  an  input  upstream  cyclic 
flow  profile  according  to  Robertson's  Recurrence  Model  for  a 
specified  range  of  alpha(a)  and  beta(/3)  values.  The  Root 
Mean  Square  Of  Error  (RMSE)  between  the  predicted  downstream 
cyclic  flow  profile  and  actual  downstream  cyclic  flow 
profile  over  all  intervals  is  calculated  to  serve  as  a 
preliminary  means  of  measuring  the  goodness-of-f i t .  This 
provides  a  relative  measure  of  the  accuracy  of  the  predicted 
cyclic  flow  profile  compared  to  the  actual  cyclic  flow 
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profile  for  a  given  street  segment.  It  should  be  noted  that 
this  type  of  comparison  would  not  be  valid  between  different 
street  segments. 

The  RMSE,  which  is  one  performance  indicator  to  measure 
the  accuracy  of  the  Robertson  model,  takes  the  following 
form  : 

n 

RMSE  =  SQRT[  I  [Actual  Flow ( i ) -Predicted  Flow(i)]2] 

i  =  1 

where  : 

RMSE  =  Root  Mean  Square  Of  Error 

i  =  1,2,...,n;  where  n=number  of  time 
intervals  comprising  the  cycle 
length  at  the  upstream  signal 

Actual  Flow  (i)  =  the  observed  flow  in  the  i-th  time 

interval  at  the  downstream  cyclic 
flow  profile,  and 

Predicted  Flow  (i)  =  the  predicted  flow  in  the  i-th  time 

interval  at  the  downstream  cyclic 
flow  profile. 

For  a  given  alpha  value,  the  performance  indicator 
RMSE,  is  calculated  for  a  complete  range  of  beta  from  a 
specified  lower  limit  to  a  specified  upper  limit  for  a  given 
increment  size.  Similarily,  for  a  given  beta  value,  the 
performance  indicator  RMSE,  is  calculated  for  a  range  of 
alpha  from  a  specified  lower  limit  to  a  specific  upper  limit 
for  a  given  increment  size.  The  analysis  for  each 
combination  of  alpha  and  beta  is  presented  in  a  tabular  and 
graphical  form.  The  combination  of  alpha  and  beta  which 
minimizes  the  RMSE  performance  indicator  is  therefore 
considered  to  be  the  'best'  calibrated  empirical  parameters. 
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A  typical  tabular  summary  of  the  G0RD0NV1  output  is 
presented  in  Figure  4.1.  The  following  information  is 
summarized : 


1.  Robertson's  beta  value,  /?T  (note  the  integer  value  in 
brackets),  alpha  value,  and  F  factor. 


2.  Four  header  cards  to  organize  the  data  base  with  respect 
to  location,  survey  data/time,  run  identifier  and 
channel  number. 


3.  A  table  with  five  columns  of  data.  Column  1  represents 
the  number  of  intervals  used  in  the  analysis  (interval 
size  times  number  of  intervals  equal  cycle  length). 
Column  2  represents  the  initial  upstream  cyclic  flow 
profile.  This  information  is  obtained  from  the  arrival 
matrix  (vehicles  per  interval  size)  for  location  #1. 
Column  3  represents  the  calculated  downstream  cyclic 
flow  profile  based  on  the  empirical  parameters  at  the 
top  of  the  table  and  Robertson's  Recurrence  Model. 
Column  4  represents  the  actual/observed  downstream 
cyclic  flow  profile.  This  information  is  obtained  from 
the  arrival  matrix  (vehicles  per  interval  size)  for 
location  #3.  Column  5  indicates  for  each  interval  the 
error  between  the  actual  downstream  and  predicted 
downstream  flows. 

In  reviewing  Figure  4.1,  several  items  are  worth  noting 
as  they  relate  to  the  data  collection  technique,  and  more 
importantly,  to  Robertson’s  Recurrence  Model. 


1.  Generally  speaking,  the  average  cyclic  volume  for  the 
initial  cyclic  flow  profile  should  be  equal  to  the 
actual  downstream  cyclic  flow  profile  if  no  overcounting 
or  undercounting  takes  place.  Using  the  information 
presented  in  Figure  4.1  we  can  tabulate  and  find  out 
that  the  average  flow  of  the  initial  cyclic  flow  profile 
is  40.97  while  the  average  flow  of  the  actual  downstram 
cyclic  flow  is  41.58.  This  discrepancy  is  due  solely  to 
an  error  in  the  data  collection  technique  (surveyor 
error ) . 


2.  Consistently,  the  average  cyclic  volume  for  the 
predicted  downstream  cyclic  flow  profile  will  be  less 
than  the  average  cyclic  volume  for  the  initial  cyclic 
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ROBERTSON'S  EOUflT  I  ON  AS  FT  EXISTS  Tn  The  TRANS  YT  MODEL 
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flow  profile.  This  error  can  be  calculated  using  the 
error  equation  presented  in  Section  2.3  where  (1-F)/F 
takes  on  the  value  ( 1 -0 . 25 ) /0 . 25  =  3.00  and 

q'  =  q' ( s ) 

(n+0T) 

takes  on  the  value  0.15.  The  product  of  these  two  values 
(3.0  x  0.15  =  0.45)  subject  to  rounding  represents  the 
absolute  error.  For  this  specific  example  we  conclude 
that  Robertson's  equation  underestimated  the  downstream 
average  cyclic  flow  profile  by  0.45  vehicles  or 
alternatively  by  a  factor  of  1.01. 


3.  The  calculated  random  delay  for  the  average  cyclic 
volume  for  the  actual  downstream  cyclic  flow  profile  is 
6.08  seconds  per  vehicle.  The  calculated  random  delay 
for  the  average  cyclic  volume  for  the  predicted 
downstream  cyclic  flow  profile  is  4.51  seconds  per 
vehicle,  or  alternately,  an  underestimation  of  the 
random  delay  by  a  factor  of  1.35.  (See  Section  2.3  and 
Sec  ton  4.3). 


4.  For  this  specific  example,  Robertson's  Recurrence  Model 
will  underestimate  the  predicted  downstream  volume  by  a 
factor  of  1.01  which  will  in  turn  result  in  an 
underestimate  of  the  random  delay  by  a  factor  of  1.35. 

Figure  4.2  illustrates  a  typical  graphical  summary  of 
the  goodness-of-f i t  between  the  predicted/actual  downstream 
cyclic  flow  profile  based  on  the  information  presented  in 
Figure  4.1.  The  following  specific  information  is  noted: 


1.  Four  header  cards  to  organize  the  data  base  with  respect 
to  location,  survey  date/time,  run  identifier  and 
channel  number. 


2.  Plot  of  the  measured/actual  downstream  cyclic  flow 
profile  which  is  indicated  by  an  asterisk  (*). 

3.  Plot  of  the  predicted  downstream  cyclic  flow  profile 
which  is  indicated  by  a  letter  (o).  A  visual  inspection 
of  these  graphs  indicate  that  during  severe  winter 
conditions  Robertson's  recommended  parameter  of 
alpha=0.5  and  beta=0.8  will  produce  significant  error. 


. 


. 
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R06fcfetS0N'S  EOUAf  ION  AS  TT  '  6  X  1  s  T  $  In  T  h  f  ft  A  N  S  V  t  MODEL - - - 

LOCATION  104  AVENUE  -  116  STREET  ;  EASTBOUND  VEHICULAR  FLOWS 

SURVEY  OATE/TIHE  :  THURSDAY  ,  FEBRUARY  4  ,  1962  ;  7:30  A.M. 

RUN  IDENTIFIER  WINTER  CONDITIONS  ,  EXHAUST  FUMES  IN  THE  AIR  ICY  PAVEMENT 

CHANNEL  NUMBER  ( O  ♦  1)  /  (  4  4  5  )  8 0  /  426  FEET  0  OWN  S  TR  E  AM  OF  STOPBAR 


MEASUREO  VS  PREDICTED  CYCLIC  FLOW  PROFILES 


-•—MEASURED  CYCLIC  FLOW  PROFILE 
0  PREDICTED  CYCLIC  FLOW  PROFILE 

ROOT  MEAN  SOUARE  OF  ERROR  IS  2.463 

■ 19:33  T*0.15S  RC«0 


FIGURE  4.2 


TYPICAL  GRAPHICAL  SUMMARY  OF  THE  GORDONV 1  OUTPUT 
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4.  The  X  or  horizontal  axis  is  the  cycle  time  in  steps.  The 
Y  or  vertical  axis  is  flow  per  unit  of  time. 

5.  The  RMSE  performance  indicator  for  the  above  two  plots 
is  summarized. 

Figure  4.1  and  4.2  are  repeated  for  all  combination  of 
alpha  and  beta  specified  by  the  user. 

Before  using  the  G0RD0NV1  program  to  calibrate 
Robertson's  Recurrence  Model,  it  is  instructive  to 
understand  the  functional  relationship  between  the  value  of 
RMSE  and  Robertson's  alpha  and  beta  parameters.  For  a  given 
beta  value,  the  range  of  RMSE  is  calculated  for  a  large 
range  of  alpha  from  0.05  to  0.60.  The  results  of  this 
analysis,  using  data  from  severe  winter  conditions,  is 

presented  in  Figure  4.3.  Several  observations  are  noted. 
Firstly  there  exists  a  unique  value  of  alpha,  which  can 
minimize  the  global  RMSE.  Secondly,  for  this  unique  value  of 
alpha,  there  are  several  values  of  beta  which  also 

contribute  to  a  global  minimal  value  of  RMSE.  One  unique 
family  of  curves  exists  for  each  set  of  yST  values  which 
correspond  to  one  integer  value.  Thirdly,  for  different 
ranges  of  beta  values,  there  also  exists  a  value  which 

minimizes  the  RMSE;  however,  this  value  represents  only  a 

local  minimum  vs  a  global  minimum.  Forthly,  an  alpha  value 
of  0.40  vs  the  standard  value  of  0.50  minimizes  the  value  of 
RMSE. 

For  a  given  alpha  value,  the  range  of  RMSE  has  also 
been  calculated  for  a  large  range  of  beta  from  0.30  to  0.80. 
The  results  of  this  analysis  are  presented  in  Figure  4.4. 


. 


ROOT  MEAN  SQUARE  OF  ERROR  (R.M.S.E. ) 
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1 1 


legend: 


2.00- 


—  BETA (  p  )  =  0.50,  0.55,0.60 

-  BETA  ( P  )  =  0.65,  0.70,  0.75 


1.80  - 


1.60  - 


DESIGN  NOTES 

AVERAGE  JOURNEY  TIME  FOR  WINTER 
DATA  EQUALS  7. 02  STEPS  ( 14.04  SECONDS  ) 

PxT= _ (INTEGER  VALUE) 

0.50x7.02=  3.511  (4) 

0.55x7.02  =  3.862  (4) 

0.60x7.02=  4.213  (4) 

0.65x7.02  =  4.564  (5) 


1.40  - 


1.20  - 


1.00  - 


0.80- 


0.10 


0.20  0.30  0.40 

ALPHA  VALUE  (^) 

FIGURE  4.3 


0.50 


0.60 


FUNCTIONAL  RELATIONSHIP  BETWEEN  THE  ALPHA  VALUE  AND 
ROOT  MEAN  SQUARE  OF  ERROR  FOR  VARIOUS  BETA  VALUES 


—i — 

0.70 


ROOT  MEAN  SQUARE  OF  ERROR  (R.M.S.E.  ) 
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FIGURE  4.4 

FUNCTIONAL  RELATIONSHIP  BETWEEN  THE  BETA  VALUE  AND 
ROOT  MEAN  SQUARE  OF  ERROR  FOR  VARIOUS  ALPHA  VALUES 
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The  following  observations  are  noted.  Firstly,  there  does 
not  exist  a  unique  value  of  beta  which  can  minimize  the 
RMSE.  This  observation  is  related  to  the  form  of  Robertson's 
equation  in  that  the  value  of  /?T  is  an  integer  value. 
Therefore,  for  the  range  of  beta  that  minimizes  the  RMSE, 
the  product  of  beta  times  travel  time  in  steps  represents  a 
unique  integer  value.  This  analysis  also  suggests  for  the 
range  of  beta  values  which  will  give  a  minimal  RMSE  error, 
there  will  also  exist  over  that  range  of  beta  values  only 
one  smoothing  factor  F.  We  conclude  that  the  selection  of 
the  best  beta  value  is  the  beta  value  which  minimizes  the 
rounding  error  of  the  product  of  beta,  times  the  travel  time 
in  seconds  (equals  an  integer  value).  Secondly,  a  beta  value 
of  0.570  (this  value  minimizes  the  rounding  error)  vs  the 
standard  value  of  0.80  minimizes  the  value  of  RMSE. 

The  effect  of  alpha  and  beta  on  RMSE  is  a  very 
important  consideration  in  the  calibration  process,  due  to 
the  complementary  relationship  between  alpha,  beta  and  the 
smoothing  factor.  This  information  highlights  some  of  the 
deficiencies  in  previous  research  efforts(see  Table  1.1).  In 
specific,  it  is  noted  that  those  studies  which  calibrated 
only  alpha  may  not  have  the  'best'  values. 

Based  on  the  above  analysis  the  following  procedure  was 
adopted  for  the  calibration  process: 

1.  Search  for  the  optimal  combination  of  alpha  and  beta 
which  will  minimize  the  RMSE. 

2.  Select  the  alpha  value  which  minimizes  the  RMSE. 


83 


3.  For  the  minimal  RMSE  select  a  beta  value  which  has  the 
least  round-off  error,  (beta  times  travel  time  in  steps 
equals  an  integer  value). 

4.  The  product  of  the  selected  alpha  and  beta  value  is 
referred  to  a  K. 

4.2  The  G0RD0NV2  Computing  Program 

The  G0RD0NV2  computing  program  is  similar  to  the 
G0RD0NV1  computing  program  with  one  exception.  The  exception 
is  that  Robertson’s  Recurrence  Model  is  modified  by  the 
introduction  of  an  initialization  equation  as  shown  in 
Section  2.3.  The  program  outputs  and  analysis  process  are 
identical  to  GORDONV1.  In  all  cases  it  should  be  noted  that 
the  best  combination  of  alpha  and  beta  with  the  VI  program 
are  identical  to  the  best  combination  with  the  V2  program. 
There  is  one  difference  in  that  V2  will  consistently  result 
in  a  lower  RMSE  when  compared  to  VI  since  V2  predicted 
downstream  volumes  are  identical  to  actual  downstream 
volumes . 

To  distinguish  the  difference  between  VI  and  V2  the 
following  comments  are  printed  on  the  top  of  each  tabular 
and  graphical  output  :  GORDONV 1  -  ROBERTSON'S  EQUATION  AS  IT 
EXISTS  IN  THE  TRANSYT  PROGRAM,  G0RD0NV2  -  MODIFIED 


ROBERTSON'S  EQUATION. 
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4.3  The  ( DSP I )  DELAY/STOP/PERFORMANCE  INDEX-OFFSET  Computing 
Program 

The  DSPI  computing  program  calculates,  for  a  given 
cyclic  flow  profile  (basic  input)  and  an  artifical  fixed 
time  traffic  signal,  several  signal  performance 
characteristics.  The  computation  of  delays,  stops,  and 
performance  index  from  a  cyclic  flow  profile  is  based  on 
relationships  developed  by  Gartner(30)  and  was  validated  on 
Bloor  Street  in  Toronto.  In  specific,  for  each  possible 
signal  offset,  the  following  performance  indicators  are 
calculated : 


1.  Average  delay  per  vehicle  (uniform  plus  random  delay 
component ) . 


2.  Average  stops  per  vehicle. 


3.  Performance  index  (weighted  combination  of  stops  and 
delay).  In  specific,  the  link  performance  index  was 
established  by  summing  the  uniform  delay  component  due 
to  the  signal  design  based  on  an  average  pattern  plus 
the  random  delay  component  plus  a  stop  penalty  factor 
multiplied  by  the  number  of  vehicle  stops.  For  this 
research  a  stop  penalty  factor  of  4  was  used  to  allow 
for  economic  loss  due  to  vehicle  stoppings  (36). 

A  typical  computing  output  for  the  first  page  is 
illustrated  in  Figure  4.5.  The  following  specific 
information  is  summarized  : 


1 .  Four  header  cards  to  organize  the  data  base  with  respect 
to  location,  survey  date/time,  run  identifier  and 
channel  number. 


2. 


Echo  of  input  data:  cycle  length,  effective  green  (green 
plus  1  second),  interval  size,  release  or  saturation 
rate  (maximum  discharge  across  a  signalized  stopline 
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l>  F  L  A  Y  i  F  F  S  E  T  ANALYSIS 


D  A  T  E  :  J  A  N  11/  1 9  f  A 


LOCATION 

SURVEY  DATF  /I  IME 
RUN  IDENTIFICATION 
CHANNEL  MIJMitP 


1 04  AVEMIE  -  1  1  6  STREET  ;  EAST!  ('ll  NO  VEHICULAR  FLOWS 
THURSDAY  /  F  c  u  r<  ||  a  R  Y  4  /  1  9>-  2  ;  7:30  A.v. 

WINTER  CONDITIONS  /  EXHAUST  F  LINES  IN  THF  AIR  /  ICY  RAVE  MEM 
CHARNEL  (4  +  5)  4  ?t  FLET  DOWNSTREAM  OF  STOPGAP  CUfSERVF.D  D  A  T  A  > 


INPUT  DATA 
CYCLF  LENGTH  =  °0  SECONDS 

EFFECTIVE  GKEFN  =  50  SECONDS 

INTERVAL  SIZE 
RELEASE  RATE 


K-F ACTOR 
ARRIVAL  DATA 


=  2  SECONDS 

=3240.  EPCU/HGUR 
=  1.80  EPCU/INTERVAl. 

=  4. 

=  0.24  0.24  C  .1  4  C .  1  9  0.24  0.43  0.48  0.43  1  .33  1  .29 

1.48  1.71  1.86  1.76  1.95  1.81  1.81  1.67  1.71  1.48 

1  .62  1  .81  1  .5  7  1  .52  1  .  71  1  .33  1  .38  1.19  1  .  24  1  .  38 

i  1  .24  0.90  C  .90  0  .29  0.14  0.24  0.2-4  0.05  0.05  0.05 

0.0  0  .05  0  .0  0  .19  0 .24 


OUTPUT 
CLMMULATIVE  ARRIVALS 
RANDOM  DEALY 

SATUR AT  I  UN  RATIO 


DATA 
=  41.58 

=  2.81  VEHICLE  HOURS  PER  HOUR 

=  6.08  SECONDS  PER  VEHICLE 

=  0.9  2-1 


FIGURE  4.5 


TYPICAL  TABULAR  SUMMARY  OF  THE  DSPI  OUTPUT  (PAGE  1/2) 
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under  the  given  traffic  environmental  conditions),  K 
factor  to  weight  the  relative  importance  of  stops/delay, 
and  arrival  data  (obtained  from  the  arrival  matrix  in 
the  G0RD0NV1  or  GORDONV2  computing  outputs). 


3.  Brief  overview  of  some  output  data  including  cumulative 
average  cyclic  arrivals,  the  TRANSYT  (V5  -  V6)  random 
delay  component  and  the  degree  of  saturation  expressed 
as  a  saturation  ratio. 

A  typical  computing  output  for  the  second  page  is 
illustrated  in  Figure  4.6.  Basic  information  extracted  from 
this  output  is  indicated  by  an  asterisk(*)  and  includes  the 
following  : 


1.  The  offset  (indicated  in  steps  vs  seconds)  for  which 
delays,  stops  or  the  performance  index  are  minimized. 


2.  The  associated  absolute  values 
performance  indicators. 

The  DSPI  computing  program  was 

the  downstream  survey  location 

characteristics  for  the  following 


of 

these  three  signal 

then 

used  tc 

d  evaluate  at 

the 

signal 

performance 

types  of 

cyclic  flow 

prof i les . 

1.  Measured/observed  data. 

2.  A  Robertson  profile  which  was  not  calibrated. 

3.  A  calibrated  Robertson  profile. 

4.  A  calibrated,  modified  Robertson  profile. 

The  principle  purpose  of  this  program  is  to  be  able  to 
quantify  the  impacts  in  signal  performance  as  a  result  of  an 
improvement  in  the  model  and/or  local  calibration.  The 
errors  in  signal  performance  by  not  calibrating  and/or 
improving  the  model  can  be  highlighted. 


"S  W  ; 
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S  II  M 

M  A  ff  Y  f 

F  A  N  A  L  Y  S 

I  S 

DATF:  JAN  11/  19«4 

INTER 

UN  I  F  f'R  M 

T  PT  A  L 

AV  DELAY  «  OF 

STOPS 

AVERAGE  ST NFS 

PERFORMANCE 

V  AL 

DTI  A Y  ( S  F  C  ) 

DELAY ( SFC  ) 

/VEHICLE  PER 

CYCLE 

/VEHICLE 

INDEX 

1 

1  no  2  . 1  5 

12  54  .f;9 

30.18 

41  .1  S 

0.99 

15.7  7 

? 

1063.81 

1316.55 

31.66 

41.10 

0.99 

16.46 
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1378.11 

33.14 

41.10 

0.99 

17.14 

A 

1195.53 
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41.20 
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17.92 
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1 261 .59 

1514.33 

36.42 

41.25 

0.99 

18.66 
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1  323  .4  5 

1576.19 

37.91 

41.15 

0.99 

19.34 
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1  368  .7  7 

1621.51 

39.00 

40.91 

0.98 

19.83 

O 

U 

1409 .3  1 

1662 . C  5 

39.97 

40.67 

0.98 

20.27 
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1453.66 

1706.41 

41  .04 

40.67 

0.98 

20.77 

1  C 

1421  .1  4 

1673.89 

40.26 

39.34 

0.95 

20 . 35 

1 1 

1391  .46 

164  4 . 21 

39.54 

38.53 

0.93 

19  .98 

1? 

1  34  5  .7  1 

1S98 .45 

38.44 

37.05 

0.89 

19.41 

13 

1280.9  1 

1533  .65 
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35.34 

0.85 

18.61 

14 

1204 . 3  5 

1457. 09 

35.04 

35.34 

0.85 

1  7.76 

1  5 

1135.43 

1  388  .  17 

33.39 

35.34 

0.85 

16.99 

1  0 

1052.67 

1305.41 

31  .40 

35.34 

0.85 

16  .08 

17 

980.31 

1233 . C5 

29.66 

35.34 

0.85 

15  .27 

18 

907.99 

1160.73 

27.92 

35.34 

0.85 

14.47 

19 

844 .9  5 

1097.69 

26.40 

35.34 

0.85 

13.77 

20 

778  .£9 

1031 .64 

24.81 

35.34 

0.85 
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FIGURE  4.6 

TYPICAL  TABULAR  SUMMARY  OF  THE  DSPI  OUTPUT  (PAGE  2/2) 
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In  summary,  the  DSPI  program  becomes  the  tool  to 
quantify  the  implications  of  calibration  or  modifications. 
In  contrast,  a  statistical  test  such  as  chi-squared,  would 
only  measure  statisticaly  how  closely  the  predicted  cyclic 
flow  profile  comes  to  the  observed  cyclic  flow  profile.  This 
latter  test,  gives  no  consideration  to  the  implications  of 
the  goodness-of - f it  from  the  perspective  of  delays,  stops 
and  performance  index. 
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5.  DISPERSION  OF  TRAFFIC  PLATOONS 


5.1  The  Process  Of  Calibrating  Robertson’s  Recurrence  Model 
The  process  of  calibrating  Robertson’s  recurrence  model 
is  presented  in  the  following  sequencial  steps: 


1.  For  the  measured/observed  downstream  cyclic  flow 
profile,  the  DSPI  computing  program  was  used  to  generate 
delay/offset,  stop/offset  and  performance  index/offset 
curves.  Based  on  these  curves,  the  following  information 
was  summarized  : 


Offset  at  which  the  actual  delay  is  a  minimum,  and 
the  true  value  of  the  delay  at  its  minimum  value. 


Offset  at  which  the  actual  stops  is  a  minimum,  and 
the  true  value  of  the  stop  at  its  minimum  value. 

Offset  at  which  the  actual  performance  index  is  a 
minimum,  and  the  true  value  of  the  performance  index 
at  its  minimum  value. 


In  reviewing  this  information,  it  should  be  recognized 
that  the  offset  value  which  minimizes  delay  does  not 
imply  that  stops  or  the  performance  index  are  also 
minimized  at  the  same  offset  value.  This  is  due  to  the 
fact  that,  essentially,  delays,  stops  and  performance 
index  are  different  measures  of  effectiveness  for  signal 
control.  This  observation  can  be  confirmed  by  referring 
to  Figure  4.6  which  indicates  that  a  different  minimium 
signal  offset  value  occurs  for  delay,  stops  and  the 
performance  index. 


2.  The  G0RD0NV1  computing  program  was  used  to  generate  two 
specific  analysis  for  each  site.  The  first  analysis 
involved  the  calculation  of  the  downstream  cyclic  flow 
profile  based  on  Robertson's  recommended  parameters  of 
alpha=0.50  and  beta=0.80  and  the  resulting  RMSE 
statistic.  The  second  analysis  involved  the  use  of  the 
G0RD0NV1  program  to  search  out  the  combination  of  alpha 
and  beta  which  reuslts  in  a  global  minimum  value  of  the 
RMSE  performance  indicator.  This  process  generated  a 
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standard  Robertson  profile  and  a  calibrated  downstream 
cyclic  flow  profile  based  on  the  Robertson's  Recurrence 
Model.  The  DSPI  computing  program  was  used  to  generate  a 
delay/stop/per formance- index  offset  infomation.  Based  on 
these  predicted  cycle  flow  profiles  and  the  actual 
cyclic  flow  profile,  the  following  statistics  were 
summarized  in  the  tabular  form  : 
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by  the  predicted  cyclic  flow  profile. 
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the  predicted  cyclic  flow  profile 
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f  the  performance 
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ted  from  the  actual 
offset  given  by  the 


For  a  perfect  fit  between  the  predicted  downstream 
cycl ic  flow  prof i le  and  the  actual  downstream  cyclic 
flow  prof i le ,  the  di f  f erences  between  the  var i ous 
performance  indicators  at  the  best  offset  value  or  over 
all  offset  values  would  be  zero. 


The  previously  described  process  was  repeated  with  the 
G0RD0NV2  computing  program  with  the  exception  that  there 
was  no  need  to  search  out  the  combination  of  alpha  and 
beta  which  results  in  a  global  minimum  value  of  the  RMSE 
performance  indicator. 
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5.2  Calibration  Results 

For  each  of  the  six  survey  sites,  four  cases  are 
analyzed  from  the  perspective  of  delays,  stops  and 
performance  index.  These  cases  are: 

Case  I  :  Robertson's  Model  (Robertson's  recommended 
alpha  value  of  0.50  and  beta  value  of  0.80  is  assumed.) 

Case  II  :  Calibrated  Robertson's  Model  (Calibration  of 
both  alpha  and  beta  values). 


Case  III  :  Calibrated,  Modified  Robertson's  Model 
(Introduction  of  an  initialization  equation  and 
calibration  of  both  alpha  and  beta  values). 

Case  IV  :  Modified  Robertson's  Model  (Robertson's 
recommended  alpha  value  of  0.50  and  beta  value  of  0.80 
i s  assumed ) . 

Tables  5.1  to  5.6  inclusive  provide  a  comparison  of 
platoon  prediction  errors  with  respect  to  delay,  stops  and 
performance  index  for  all  intersection  surveyed.  To  explain 
in  more  detail  the  exact  process  and  the  results  obtained, 
Table  5.1  will  be  explained  in  detail. 

For  the  measured/observed  downstream  cyclic  flow 
profile  the  minimum  delay  occurs  at  an  offset  of  68  seconds 
and  the  true  value  of  this  delay  at  its  minimum  value  is 
9.70  seconds  per  vehicle.  Similarily,  the  minimum  stops 
occurs  at  an  offset  of  56  seconds  and  the  true  value  of  the 
stops  at  its  miminal  value  is  0.64  stops  per  vehicle.  Also, 
the  minimum  performance  index  occurs  at  an  offset  of  64 
seconds  and  the  true  value  of  the  performance  index  at  its 
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minimum  value  is  5.88. 

For  the  predicted  downstream  cyclic  flow  profile  which 
is  based  on  Robertson ’ s  Model  (alpha=0.50,  beta=0.80),  the 
minimum  delay  occurs  at  an  offset  of  72  seconds  and  the 
predicted  value  of  the  delay  at  its  minimum  value  is  8.52 
seconds  per  vehicles.  Using  the  observed  downstream  cyclic 
flow  profile,  the  true  delay  at  the  72  second  offset  value 
is  12.58  seconds  per  vehicle.  Similarily,  the  minimum  stops 
occurs  at  an  offset  of  64  seconds  and  the  predicted  value  of 
the  stops  at  its  minimum  value  is  0.30  stops  per  vehicle 
Using  the  observed  downstream  cyclic  flow  profile  the  true 
value  of  stops  at  the  64  second  offset  value  is  0.74.  Also, 
the  minimal  performance  index  occurs  at  an  offset  of  68 
seconds  and  the  predicted  value  of  the  performance  index  at 
its  minimal  value  is  4.94.  Using  the  observed  downstream 
cyclic  flow  profile,  the  true  performance  index  at  the  68 
second  offset  value  is  6.08.  In  summary,  for  this  particular 
site  with  the  use  of  Robertson’s  recommended  empirical 
parameters  (alpha=0.50,  beta=0.80)  the  delay  indicator  will 
be  in  error  by  48%,  the  stops  indicator  will  be  in  error  by 
147%  and  the  performance  indicator  will  be  in  error  by  23%. 
One  therefore  concludes  that  significant  errors  in  the 
signal  performance  characteristics  can  materialize  by  not 
calibrating  Robertson's  Recurrence  Model. 

For  the  predicted  downstream  cyclic  flow  profile  which 
is  based  on  a  Calibrated  Robertson ' s  Model  (alpha=0.40, 


beta=0.57)  we  note  that  RMSE  index  is  decreased  from  2.463 


. 
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to  0.863,  suggesting  that  a  better  fit  is  being  obtained  by 
calibration.  A  review  of  the  three  basic  signal  performance 
indicators  for  this  case  shows  that  with  calibration  the 
delay  indicator  will  be  in  error  by  24%;  however,  the 
improvement  in  the  estimate  of  delay  is  in  the  order  of  50% 
(( 48-24 ) /48 ) .  Similarily  for  stops,  the  stop  indicator  will 
be  in  error  by  31%;  however,  the  improvement  in  the  estimate 
of  stops  is  in  the  order  of  79%  ( ( 1 47-3 1 ) /I  47 ) .  Also,  for 
the  performance  index,  the  indicator  will  be  in  error  by 
24%,  however  this  makes  no  significant  change  (slight 
decrease  in  performance)  as  a  result  of  calibration. 

For  the  predicted  downstream  cyclic  flow  profile  which 
is  based  on  a  Calibrated  Modi f i ed  Robertson ’  s  Model 
(alpha=0.40,  beta=0.57),  we  note  a  further  improvement  in 
the  RMSE  index  from  0.863  to  0.844,  suggesting  a  even  better 
fit  can  be  obtained  by  providing  an  initialization  equation. 
A  review  of  the  three  basic  signal  performance  indicators 
for  this  case  shows  that  with  modifications  and  calibration 
the  delay  indicator  will  still  be  in  error  by  16%,  however 
this  further  improvement  (calibrated  modified  Robertson’s 
model  in  comparsion  to  a  calibrated  Robertson's  model)  in 
the  estimate  of  delay  is  in  the  order  of  33%  ( (24— 16)/24) . 
Similarily  for  stops,  the  stop  indicator  will  still  be  in 
error  by  28%;  however,  this  further  improvement  in  the 
estimate  of  stops  is  in  the  order  of  10%  ( (31—28) /3 1 ) .  Also 
for  the  performance  index,  the  performance  indicator  will 
still  be  in  error  by  16%;  however,  this  further  improvement 
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in  the  estimate  of  the  performance  index  is  in  the  order  of 
33%  ( ( 24- 1 6 )/3 1 ) . 

The  last  set  of  analysis  is  based  on  the  predicted 
cyclic  flow  profile,  which  is  based  on  an  Modified 
Robertson ' s  Model  (alpha=0.50,  beta=0.80).  Generally 

speaking,  in  comparison  to  the  Robertson's  Model  there  is  a 
consistent  and  significant  improvement  with  respect  to 
delays,  stops,  performance  index  and  the  RMSE  statistics. 
However  the  improvements  in  performance  are  not  comparable 
to  those  obtained  by  calibrating  and  modifing  Robertson's 
model . 

To  obtain  a  comprehensive  overview  of  all  the  analysis 
that  was  undertaken,  a  linear  regression  analysis  of  all 
predicted  signal  performance  characteristics  in  comparison 
to  the  actual  performance  characteristics  was  undertaken  in 
order  to  deduce,  from  a  statistical  perspective,  the 
improvements  relating  to  calibration  and/or  modifications  to 
the  model. 

5.2.1  Delay  Considerations 

A  regression  analysis  of  the  predicted  minimal  delay  vs 
the  actual  minimal  delay  was  undertaken  for  the  data  that 
was  analyzed  at  all  survey  sites.  The  linear  regression 
equation  was  of  the  form: 

Yd  =  ai  +  b,Xd 

where  Yd  =  predicted  minimal  delay  (seconds 

per  vehicle ) , 
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Xd  =  actual  minimal  delay  (seconds  per 
vehicle ) ,  and 

airb1  =  calibrated  coefficents. 

For  a  perfect  fit  between  predicted  minimum  delay  vs 
actual  minimum  delay,  the  coefficient  a!  would  assume  a  0.00 
value  and  the  coefficient  b,  would  assume  a  1.00  value. 
Alternately,  a  good  platoon  dispersion  model  should  approach 
these  specific  values  for  the  coefficients. 

The  resulting  regression  equation  for  CASE  I,  CASE  II 
and  CASE  III  along  with  R-Square  (R2),  and  correlation 
coefficient  (p)  are  listed  below  : 

CASE  I 

Robertson's  Model 

Y d  =  0.607  +  0 . 7 9 1 Xd  ;  R2  =  0.92,  p  =  0.96 

CASE  II 

Calibrated  Robertson's  Model 

Y d  =  0.497  +  0 . 828Xd  ;  R2  =  0.99,  p  =  0.99 

CASE  III 

Calibrated  Modified  Robertson's  Model 

Y d  =  0  .  106  +  0 . 949Xd  ;  R2  =  0.99,  p  =  0.99 

The  above  equations  indicate  that,  from  a  statistical 
perspective,  the  Calibrated  Modified  Robertson  Model  results 
in  the  best  overall  form  of  model  for  predicting  minimal 
vehicular  delay.  Figure  5.1  illustrates  all  three  delay 
regression  lines  along  with  the  data  points  for  CASE  III.  In 
reviewing  this  figure,  it  should  be  noted  that  all  five  of 
the  six  data  points  fall  within  the  95%  confidence  band.  One 
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legend: 

■  -  DATA  POINTS  FOR  THE  CALIBRATED  MODIFIED 

ROBERTSON  MODEL. 

-  -  REGRESSION  LINE  FOR  THE  CALIBRATED 

MODIFIED  ROBERTSON  MODEL. 

.  -  REGRESSION  LINE  FOR  THE  CALIBRATED 

ROBERTSON  MODEL. 

- -  REGRESSION  LINE  FOR  ROBERTSON 

MODEL.  (Not  Calibrated) 


4  5  6  7  8  9  10  II  12  13  14 

ACTUAL  MINIMUM  DELAY  (  Seconds  per  Vehicle) 


FIGURE  5. 1 

GOODNESS-OF-FIT  FOR  THREE  FORMS  OF  ROBERTSON'S  RECURRENCE 
MODEL  FROM  A  MINIMUM  DELAY  PERSPECTIVE 
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data  point  misses  this  confidence  band  by  approximately  0.2 
seconds  per  vehicle. 

For  comparison  purposes  only,  Lam's  data  l22)  was 
subjected  to  similar  regression  analysis.  The  equations 
derived  for  CASE  I  and  CASE  III  respectively  were  : 

CASE  I 

Robertson's  Model 

Yd  =  1.091  +  0 . 79 1 Xd  ;  R2  =  0.69,  p  =  0.83 

CASE  III 

Calibrated  Modified  Robertson's  Model 
Y d  =  0.997  +  0 . 9  1 2Xd  ;  R2  =  0.79,  p  =  0.89 


Principal  conclusions  reached  from  a  mimimum  delay 
perspective  are  as  follows  : 


1.  In  the  range  of  0-5  seconds  per  vehicle  the  effect  of 
calibration  or  modifications  to  Robertson's  Recurrence 
Model  is  insignificant. 


2.  In  the  range  of  5  -  15  seconds  per  vehicle,  a 

Robertson's  Model  which  has  not  been  calibrated  could  be 
in  error  up  to  20%  ((15.00/12.47)). 


3.  In  the  range  of  5-15  seconds  per  vehicle  a  Calibrated 
Robertson's  Model  could  be  in  error  up  to  16% 
( ( 15.00/12.92) ) . 


4.  In  the  range  of  5-15  seconds  per  vehicle  a  Calibrated 
Modified  Robertson's  Model  could  still  be  in  error  up  to 
7%  ((15.00/14.03)). 


5.  From  the  data  collected  and  analyzed,  no  conclusion  can 
be  reached  for  delays  greater  than  15  seconds  per 
vehicle.  However,  given  the  R2  and  p  values,  we  can 
conclude  that  this  equation  is  strong  in  its  predictive 
power  and  that  similar  conclusion  can  be  drawn  for 
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larger  delay  values. 

Figure  5.2  illustrates  graphically  the  effect  of 
calibration  or  modifications  from  a  minimum  delay 
perspective . 

5.2.2  Stop  Considerations 

Undertaking  a  similar  regression  anlaysis  for  stops  as 
was  presented  in  Section  5.2.1,  the  following  regression 
equations  were  obtained: 

CASE  I 

Robertson's  Model 

Y  s  =  0  .  1  05  +  0 . 4  3  3X  s  ;  R2  =  0.74,  p  =  0.74 

CASE  II 

Calibrated  Robertson's  Model 

Y  s  =  0.000  +  0 . 702X  s  ;  R2  =  0.95,  p  =  0.97 

CASE  III 

Calibrated  Modified  Robertson's  Model 

Y s  =  0.000  +  0 . 7 2  1 X $  ;  R2  =  0.94,  p  =  0.96 

Based  on  the  above  equations,  we  conclude  that,  from  a 
statistical  perspective  the  Calibrated  Modified  Robertson 
Model  results  in  the  best  overall  form  of  the  model  for 
predicting  minimum  vehicular  stops.  Figure  5.3  illustrates 
all  three  stop  regression  lines  along  with'  the  data  points 
for  CASE  III.  In  reviewing  this  figure  it  should  be  noted 
that  all  six  data  points  fall  within  the  95%  confidence 


bands . 


' 


RELATIVE  INDEX 


105 


FIGURE  5.2 
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EFFECT  OF  CALIBRATION  OR  MODIFICATIONS  FROM  A  MINIMUM  DELAY 

PERSPECTIVE 


ERROR  IN  THE  PREDICTION 


PREDICTED  MINIMAL  STOPS  (Stops  per  Vehicle  ) 
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legend: 
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MODIFIED  ROBERTSON  MODEL. 
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FIGURE  5.3 

GOODNESS-OF-FIT  FOR  THREE  FORMS  OF  ROBERTSON'S  RECURRENCE 
MODEL  FROM  A  MINIMUM  STOPS  PERSPECTIVE 
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Principal  conclusions  that  are  reached  from  a  minimum 
stop  perspective  are  as  follows  : 


In  the  range  of  0.35-0.45  stops  per  vehicle  the  effect 
of  calibration  of  modifications  to  Robertson's  model  is 
ins ign i f icant . 


2.  In  the  range  of  0.45-0.80  stops  per  vehicle  a 
Robertson's  Model  which  has  not  been  calibrated  could  be 
in  error  by  78%  ( ( 0 . 80 )/0 . 45) ) .  In  the  range  of 
0.00-0.45  stops  per  vehicle,  a  Robertson's  Model  which 
has  not  been  calibrated  could  have  larger  percentage 
errors;  however,  the  magnitude  of  this  error  in  absolute 
values  is  smaller.  The  former  errors  represent  an 
underestimation  while  the  later  errors  represent  an 
overest imat i on . 


3.  A  Calibrated  Robertson's  Model  could  be  in  error  up  to 
43%  ((0.80/0.56)). 


4.  A  Calibrated  Modified  Robertson's  Model  could  be  in 
error  up  to  38%  ( ( 0 . 80 )/0 . 58 ) ) . 

5.  From  the  data  collected  and  analyzed,  no  conclusions  can 
be  reached  for  minimum  stops  greater  than  0.85  stops  per 
vehicle.  However,  given  the  R2  and  p  values,  we  can 
conclude  that  this  equation  is  strong  in  terms  of  its 
predictive  power  and  that  similar  conclusions  can  be 
drawn  for  larger  stop  values. 

Figure  5.4  illustrates  graphically  the  effect  of 
calibration  and  modifications  from  a  minimum  stop 
perspective . 


5.2.3  Performance  Index  Considerations 


Again,  a  similar  regression  analysis  for  the 
performance  index  as  was  presented  in  Section  5.2.1  and 
Section  5.2.2  resulted  in  the  following  regression 
equations: 
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FIGURE  5.4 

EFFECT  OF  CALIBRATION  OR  MODIFICATIONS  FROM  A  MINIMUM  STOPS 

PERSPECTIVE 


ERRORS  IN  THE  PREDICTIONS  IN  COMPARISON  TO  THE  MEASURED/ 
OBSERVED  DATA. 
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CASE  I 

Robertson's  Model 

Y  =  0.041  +  0.877X  ;  R2  =  0.99,  p  =  0.99 

P  P 

CASE  II 

Calibrated  Robertson's  Model 

Y  =  0.120  +  0 . 86 IX  ;  R2  =  0.99,  p  =  0.99 

P  P 

CASE  III 

Calibrated  Modified  Robertson's  Model 

Y  =  0.003  +  0.963X  ;  R2  =  0.99,  p  =  0.99 

P  P 

Based  on  the  above  equations,  we  can  conclude,  that 
from  a  statistical  perspective,  the  Calibrated  Modified 
Robertson  Model  results  in  the  overall  best  form  of  the 
model  for  predicting  minimum  performance  index.  Figure  5.5 
illustrates  all  three  performance  index  regression  lines 
along  with  the  data  points  for  CASE  III.  In  reviewing  this 
figure,  it  should  be  noted  that  all  five  of  the  six  data 
points  fall  within  the  95%  confidence  band.  One  data  point 
misses  this  confidence  interval  by  only  0.09. 

Principle  conclusions  that  can  be  reached  from  a 
minimum  performance  index  perspective  are  as  follows  : 

1.  Robertson's  Model  which  has  not  been  calibrated  could  be 
in  error  up  to  14%  ( ( 1 0 . 00 ) /8 . 82 ) ) . 

2.  A  Calibrated  Robertson's  Model  could  be  in  error  up  to 
15%  ( ( 10.00) /8 .73).  The  differences  between  a  calibrated 
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legend: 


■  -  DATA  POINTS  FOR  THE  CALIBRATED  MODIFIED 

ROBERTSON  MODEL. 

—  -  REGRESSION  LINE  FOR  THE  CALIBRATED 


ACTUAL  PERFORMANCE  INDEX 


FIGURE  5.5 

GOODNESS-OF-FIT  FOR  THREE  FORMS  OF  ROBERTSON'S  RECURRENCE 
MODEL  FROM  A  MINIMUM  PERFORMANCE  INDEX  PERSPECTIVE 


and  a  non-cal ibrated  Robertson's  Model  appear  to  be 
insigni f icant . 

3.  A  Calibrated  Modified  Roberton's  Model  could  be  in  error 
up  to  by  4%  ( (  1  0 .00  )/96 .4  )  . 

Figure  5.6  illustrates  graphically  the  effect  of 
calibration  and  modifications  from  a  minimum  performance 
index  perspective.  It  should  be  noted  that  the  use  of  a 
different  weighting  factor  (performance  index  equals  total 
delay  (uniform  plus  random)  plus  four  times  the  number  of 
stops)  may  alter  these  conclusions. 

5.2.4  Overall  Impact  Of  Calibration  And  Modifications  To 
Robertson's  Model 

Based  on  the  previously  discussed  regression  analysis 
we  conclude  that  the  Calibrated  Modi f ied  Robertson ' s  Model 
can  improve  the  predictive  estimates  of  the  minimum  value  of 
delays  by  13%  (93.5/83.1),  stops  by  29%  (72.5/56.3),  and  the 
performance  index  by  9%  (96.4/88.2).  On  the  other  hand,  a 
Calibrated  Robertson ' s  Model  can  improve  on  the  previously 
defined  signal  performance  indicators  by  4%,  24%  and  11% 
respectively.  A  similar  regression  analysis  was  also 
duplicated  for  the  complete  range  of  offsets  (all  offset 
values  including  the  minimum  value )  for  delays  and  stops  at 
each  intersection  surveyed  and  analyzed.  Similar  conclusive 
results  were  obtained.  Calibration  of  the  individual  value 
of  alpha  and  beta  along  with  some  modifications  to  the 
model,  will  consistently  improve  the  prediction  of  signal 
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FIGURE  5.6 

EFFECT  OF  CALIBRATION  OR  MODIFICATIONS  FROM  A  MINIMUM 

PERFORMANCE  INDEX  PERSPECTIVE 


OBSERVED  DATA. 
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performance  indicators  such  as  delays,  stops  and  performance 
index  for  any  signal  offset  value. 


5.3  Design  Guidelines 

The  calibration  results  for  the  survey  sites 
investigated  are  summarized  in  Table  5.7,  and  illustrated 
graphically  in  Figures  5.7  to  5.12.  For  each  site 
investigated  Table  5.7  presents  the  following  information  : 
the  alpha  value  which  minimized  the  RMSE,  the  beta  value 
which  minimized  both  the  RMSE  and  the  rounding  off  error 
when  the  product  of  beta  and  the  travel  time  in  steps  is 
calculated,  the  resulting  K  value  (note  Robertson's 
recommended  K  value  is  40),  the  calculated  smoothing  factor 
F,  and  the  associated  degree  of  saturation.  In  reviewing 
this  table  the  following  trends  have  been  identified  : 


ALPHA  VALUE 


1.  At  intersections  (summer  data)  where  the  degree  of 
saturation  ranges  from  a  high  value  (104  Avenue  -  116 
Street,  Week  No.  1  and  2)  to  a  medium  value  (111  Avenue 
-  116  Street)  to  a  low  value  (102  A  Avenue  -  97  Street), 
we  note  a  corresponding  decrease  in  the  value  of  alpha 
ranging  from  0.45  to  0.15 


2.  The  basic  difference  between  summer  and  winter 

conditions  at  high  degrees  of  saturation  appear  to  be  a 
slight  decrease  in  the  alpha  value  (Summer  0.45  to  0.40 
vs  Winter  0.40  to  0.30).  Severe  winter  conditions  are  at 
the  higher  end  of  the  range  while  average  winter 
conditions  are  at  the  lower  end  of  the  range. 


CALIBRATED  PARAMETERS  FOR  THE  SURVEY  SITES  INVESTIGATED  DURING  THE  AM  PEAK  HOUR 
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4  .  S  I 

4.4  I 

4.3  I 

4.2  1 

4.1  I 

4.0  I 
3  9  I 

3  8  I 

3.7  I 

3.6  I 

3  .  5  1 

3.4  I 

3.3  I 

3.2  I 

3  1  I 

3  0  I 

2  9  I 

2.6  I 

2.7  I 
2.6  I 
2.6  I 

2.4  I 

2.3  I 
2.2  I 
2.1  I 
2.0  I 


— ^-MEASURED  CYCLIC  FLOW  PROFILE 
0  PREDICTED  CYCLIC  FLOW  PROFILE 


pOOT  MEAN  SQUARE  OF  ERROR  IS_ 0.644 


FIGURE  5.7 

COMPARISON  OF  MEASURED  (*)  VS  PREDICTED  (0)  CYCLIC  FLOW 
PROFILE  FOR  104  AVENUE  -  116  STREET  (SEVERE  WINTER 
CONDITIONS)  WITH  THE  FOLLOWING  CALIBRATED  ROBERTSON’S 
PARAMETERS  :  ALPHA  =  0.40,  BETA  =  0.570 
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FIGURE  5.8 

COMPARISON  OF  MEASURED  (*)  VS  PREDICTED  (0)  CYCLIC  FLOW 
PROFILE  FOR  104  AVENUE  -  116  STREET  (AVERAGE  WINTER 
CONDITIONS,  FURTHER  DOWNSTREAM)  WITH  THE  FOLLOWING 
CALIBRATED  ROBERTSON'S  PARAMETERS  :  ALPHA  =  0.30,  BETA  = 
0.725 


MODIFIED  ROBERTSON'S  EOuAt  1  ON 


|L  OC  AT  I  ON 

(survey  DATE/TIME 
(run  IDENTIFIER 
CHANNEL.  NUMBER 


104  AVENUE  -  IIS  STREET  ;  EASTBOUND  VEHICULAR  FLOWS 
THURSDAY  APRIL  16  ,  1981  ;  7:30  A.M. 

SUMMER  CONDITIONS  ,  DRY  PAVEMENT 

<0  ♦  1)  /  (4  ♦  S)  ;  80  /  426  FEET  DOWNSTREAM 


4  .  S 
4  .  4 
4  .  3 
4  .  2 
4  .  1 
4  .  O 
3 . 9 
3  .  8 
3  .  7 
3  .  6 
3  .  5 
3  .  4 
3  .  3 
3 . 2 
3  .  1 
3  .  O 
2  .  9 
2 . 8 
2 . 7 
2  .  9 


MEASURED  VS  PREDICTED  CYCLIC  FLOW  PROFILES 


CYCLE  TIME  (  2  SECOND  INTERVALS)  FROM  START  OF  MAIN  STREET  GREEN 

---MEASURED  CYCLIC  FLOW  PROFILE 
0  'PREDICTED  CYCLIC  FLOW  PROFILE 


ROOT  MEAN  SOUARE  OF  ERROR  IS 


1.009 


FIGURE  5.9 

COMPARISON  OF  MEASURED  (*)  VS  PREDICTED  (0)  CYCLIC  FLOW 
PROFILE  FOR  104  AVENUE  -  116  STREET  (WEEK  NO.  1)  WITH  THE 
FOLLOWING  CALIBRATED  ROBERTSON'S  PARAMETERS  :  ALPHA  =  0.45, 
BETA  =  0.811 
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MEASURED  VS  PREDICTED  CYCLIC  FLOW  PROFILES 
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ROOT  MEAN  SQUARE  OF  ERROR  IS 
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FIGURE  5.10 

COMPARISON  OF  MEASURED  (*)  VS  PREDICTED  (0)  CYCLIC  FLOW 
PROFILE  FOR  104  AVENUE  -  116  STREET  (WEEK  NO.  2)  WITH  THE 
FOLLOWING  CALIBRATED  ROBERTSON'S  PARAMETERS  :  ALPHA  =  0.40, 
BETA  =  0.811 
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MEASURED  VS  PREDICTED  CYCLIC  FLOW  PROFILES 


CYCLE  TIME  (  2  SECOND  INTERVALS)  FROM  START  OF  MAIN  STREET  GREEN 

-•-MEASURED  CYCLIC  FLOW  PROFILE 
0  PREDICTED  CYCLIC  FLOW  PROFILE 


ROOT  MEAN  SQUARE  OF  ERROR  IS 


0.640 


FIGURE  5.11 

COMPARISON  OF  MEASURED  (*)  VS  PREDICTED  (0)  CYCLIC  FLOW 
PROFILE  FOR  111  AVENUE  -  116  STREET  WITH  THE  FOLLOWING 

CALIBRATED  ROBERTSON'S  PARAMETERS  :  ALPHA  =  0.35,  BETA  = 
0.722 
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-•-MEASURED  CYCLIC  FLOW  PROFILE 
0  PREDICTED  CYCLIC  FLOW  PROFILE 
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FIGURE  5.12 

COMPARISON  OF  MEASURED  (*)  VS  PREDICTED  (0)  CYCLIC  FLOW 
PROFILE  FOR  102  A  AVENUE  -  97  STREET  WITH  THE  FOLLOWING 
CALIBRATED  ROBERTSON'S  PARAMETERS  :  ALPHA  =  0.15,  BETA  = 
0.627 
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BETA  VALUE 


1.  At  intersections  (summer  data)  where  the  degree  of 
saturation  ranges  from  a  high  value  (104  Avenue  -  116 
Street  Week  No.  1  and  2)  to  a  medium  value  (111  Avenue  - 
116  Street)  to  a  low  value  (102  A  Avenue  -  97  Street), 
we  note  a  corresponding  decrease  in  the  value  of  beta 
ranging  from  0.811  to  0.627. 


2.  The  basic  difference  between  summer  and  winter 
conditions  at  high  degrees  of  saturation  appears  to  have 
significant  impact  on  the  range  of  beta  values.  During 
typical  winter  conditions,  beta  is  slightly  lower  than 
summer  conditions;  however,  during  severe  winter 
conditions,  the  beta  value  observed  is  lower  than  the 
lowest  value  of  all  beta  values  in  summer. 

K  VALUE 


1.  When  the  TRANSYT  model  was  initially  introduced,  the 
recommended  K  value  was  40  (alpha=0.50,  beta=0.80).  The 
most  current  version  of  TRANSYT,  Version  8,  recommends  a 
K  value  of  35  (alpha=0.44,  beta=0.80)  in  response  to  the 
available  published  literature  which  indicates  that  more 
dispersion  is  actually  taking  place  in  comparsion  to 
Robertson's  orignal  research.  Based  on  the  Edmonton  data 
a  K  value  of  35  still  appears  to  be  too  high,  suggesting 
that  in  comparsion  to  other  traffic  flow  patterns, 
Edmonton  drivers  do  not  exhibit  a  compact  platoon 
behaviour.  Edmonton  data  results  in  a  K  value  ranging 
from  9  to  37. 

Robertson  (13)stated  that  it  would  be  reasonable  to 
expect  that  the  value  of  F  be  a  function  of  various  site 
factors  such  as  road  width,  gradient,  parking  activity, 
opposing  vehiclar  flow  level,  traffic  composition,  etc.,  and 
not  simply  the  average  journey  time  (T) .  Robertson  when 


making  this  statement, 
value  would  only  be  a 
since  the  basic  input 


was  in  fact  indicating  that  the  alpha 
function  of  these  characteristics 
into  the  TRANSYT  model  (Card  Type  31  - 


Link/Node  Relation)  permits  the  user  to  input  K.  A  closer 


examination  indicates  that  through  K,  the  user  could  only 
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change  alpha  and  not  the  value  of  beta.  The  value  of  beta  = 
0.80  could  not  be  changed  by  the  user. 

This  Researcher,  with  the  same  reasoning,  suggests  that 
it  would  also  seem  reasonable  to  expect  that  /?T#  which  in 
effect  is  the  average  arrival  time  of  leading  vehicles  in 
the  platoon,  should  depend  to  some  extent  on  these  same  site 
factors  (perhaps  more  importantly  to  ambient  weather 
conditions)  and  not  simply  0.8  of  the  average  journey  time 
(T) .  The  Reseracher  further  suggests,  that  the  amount  of 
platoon  dispersion  (which  is  a  function  of  alpha,  beta  and 
travel  time)  could  be  correlated  with  the  degree  of 
saturation  at  the  upstream  signalized  intersection.  As 
degree  of  saturation  increases,  the  freedom  to  maneouvre 
decreases,  and  hence  we  would  expect  that  the  amount  of 
dispersion  to  decrease. 

The  Edmonton  data  as  presented  in  Table  5.7  provides 
some  conformation  that  alpha  and  beta  could  be  a  function  of 
site  characteristics.  It  appears  that  degrees  of  saturation 
(ratio  of  volume  to  capacity)  are  perhaps  a  better  indicator 
of  the  most  appropriate  value  of  alpha  and  beta  for 
design/analysis  purposes.  For  this  reason,  it  is  recommended 
that  the  information  presented  in  Table  5.7  be  used  as  a 
design  guideline.  The  Edmonton  data  also  suggests  that  when 
winter  signal  timing  plans  are  being  prepared,  estimates  of 
winter  saturation  flow  rates  ( 3 5 )  and  winter  travel  times 
alone,  are  not  adequate  inputs  for  signal  design.  Different 
emperical  values  of  alpha  and  beta  in  comparsion  to  summer 


' 
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conditions  are  also  required.  Based  on 
(weighted  in  proportion  to  the  number  o 
collected)  the  following  average  empir 
also  be  used  for  general  design  purposes 


a  weighted  average 
f  cycles  of  data 
ical  parameters  may 


1.  AVERAGE  EMPIRICAL  PARAMETERS  FOR  USE  IN  ROBERTSON'S 
RECURRENCE  MODEL  DURING  SUMMER  CONDITIONS  IN  EDMONTON, 
ALBERTA,  CANADA  :  ALPHA  =  0.34,  BETA  =  0.75. 


2. 


AVERAGE  EMPIRICAL  PARAMETERS  FOR  USE  IN  ROBERTSON'S 
RECURRENCE  MODEL  DURING  WINTER  CONDITIONS  I  EDMONTON, 
ALBERTA,  CANADA  :  ALPHA  =  0.34,  BETA  =  0.66. 


5.4  Practical  Applications 

Knowledge  of  platoon  dispersion  has  many 
appl icat ions ( 2 6 , 2 2 )  apart  from  its  main  use  in  deriving 
accurate  downstream  arrival  patterns  for  signal  offset 
design  and  evaluation  through  the  computation  of  delay, 
stops,  and  performance  index  relationships.  The  height  of 
the  peaks  and  trough  of  the  cyclic  flow  profile,  and  their 
shape  and  location  in  time,  provide  the  transportation 
engineer  with  a  great  deal  of  addition  information  as 
outlined  below: 


1.  Average  flow.  The  area  of  the  cyclic  flow  profile  is 
proportional  to  the  average  flow.  When  inspecting 
several  such  graphs,  sites  with  heavy  flows  are  easily 
ident i f ied . 


2.  Major/mi  nor  road  inflows.  If  the  green  times  along  with 
the  signal  phasing  sequence  at  the  upstream  intersection 
are  known,  then  the  major  and  minor  inflow  rates  can  be 
checked  by  observing  the  average  flow  during  the 
appropriate  part  of  the  cycle. 

3.  Benefits  of  signal  co-ordination  The  variation  in 


the 
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height  of  the  cyclic  flow  profile  is  a  measure  of  the 
need  for  signal  co-ordination.  If  the  height  is  nearly 
constant  during  the  cycle,  then  there  is  no  benefit  to 
be  gained  by  co-ordinating  the  upstream  and  downstream 
signals.  This  is  true  even  when  the  average  flow  is 
large.  Conversely,  a  cyclic  flow  profile  with  a 
pronounced  peak  corresponding  to  a  well  defined  platoon 
presents  a  strong  potential  for  good  signal 
co-ordination . 


4.  Saturation  flow.  Saturation  flow  is  defined  as  the  rate 
of  flow  across  the  stopline  from  a  queue,  but  it  can  be 
deduced  from  measurements  at  a  site  just  downstream  of 
an  intersection.  During  some  part  of  the  signal  cycle 
(during  that  portion  of  the  cycle  when  a  stationary 
queue  exists),  traffic  will  discharge  at  the  saturation 
rate.  With  knowledge  of  upstream  signal  timings,  the 
magnitude,  duration,  and  effect  of  length  of  green  on 
saturation  flow  can  be  checked  from  the  cyclic  flow 
prof i le . 


5.  Spare  capacity.  The  part  of  the  entire  green  time  during 
which  the  actual  flow  is  below  saturation  flow  levels 
represents  the  space  capacity  at  the  upstream  signal 
approach . 


6.  Platoon  speed  and  dispersion.  The  location  in  time  of 
the  peaks  and  troughs  of  the  cyclic  flow  profile  will 
depend  on  the  upstream  signal  timing,  the  time  of 
traffic  platoons  arriving  at  the  upstream  signal,  and 
the  distance  between  the  up-stream  signals  and  the  site 
at  which  the  cyclic  flow  profile  is  measured. 


7.  Signal  start  and  end  lags.  During  a  fully  saturated 
signal  phase  the  observed  cyclic  flow  profile  can  be 

used  to  calculate  signal  start  and  end  legs  (31).  A 

start  lag  is  the  initial  hesitation  following  the 
display  of  green  prior  to  traffic  discharging  at  a 

constant  rate.  A  end  lag  is  the  time  when  the  display  of 

amber  results  in  traffic  discharging  at  a  rate  which 
reduces  sharply  to  zero. 


8.  Incident  detection.  The  shape  of  the  cyclic  flow  profile 
will  alter  during  the  normal  day.  It  may  alter  more 
rapidly  if  there  is  an  accident  or  other  incident  at  an 
upstream  signal;  change  in  cyclic  flow  profile  shape  may 
be  a  sensitive  indicator  of  traffic  incidents.  (This 
specific  type  of  research  is  presently  being  undertaken 
jointly  by  the  Civil  Engineering  and  Electrical 
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Engineering  Departments  at  the  University  Of 
Alberta . ) ( 2 8 )  The  computing  program  BRAD#+*PLOTLI B 
permits  the  user  to  view  graphically  the  variations  (see 
Figure  3.4)  in  the  cyclic  flow  profile  from  one  cycle  to 
the  next. 


9.  Signal  setting  by  TRANSYT .  Cyclic  flow  profile  traffic 
analysis  is  likely  to  be  particularly  useful  when  the 
TRANSYT  method  of  optimizing  fixed-time  plant  has  been 
employed  to  derive  the  signal  timings.  The  TRANSYT 
computer  program  prints  a  cyclic  flow  profile  graph  for 
each  part  of  the  network.  Comparision  of  predicted  and 
measured  cyclic  flow  profile's  can  be  inspected  and 
quickly  reveal  any  serious  differences  between  theory 
and  practice.  Corrections  can  then  be  made  to  the  data 
fed  to  TRANSYT  and  an  improved  standard  of  traffic 
control  will  result. 


10.  International  exchange  standard.  in  the  view  of 

Robertson,  there  is  considerable  merit  in  using  cyclic 
flow  profiles  to  exchange  information  on  the  behaviour 
of  urban  road  traffic  in  different  cities  and  countries. 
Comparison  between  measured  and  predicted  cyclic  flow 
profiles  confirms  the  accuracy  of  the  existing  methods 
and  highlights  areas  for  further  research  and 

development.  For  example,  advisory  speed  are  displayed 
to  motorists  in  some  co-ordinated  signal  systems.  To  be 
effective  such  signs  would  be  expected  to  alter  the 
downstream  CFP  shape.  It  would  be  interesting  to 
establish  whether  they  do. 

As  a  further  example,  it  may  be  possible  as  this 
reseach  has  indicated  on  a  preliminary  level  to 
establish  relationships  between  weather  conditions  and 
the  cyclic  flow  profile  shape.  An  icy  road  not  only 
affects  saturation  flow  rates,  but  also  journey  times 
and  platoon  dispersion  factors.  At  the  present  time, 
little  is  known  of  the  magnitudes  of  these  effects. 

Teply  (32)  provided  a  definition  of  a  framework  for 
evaluation  of  the  applicability  of  foreign  experience  in 
traffic  management,  considering  behavioural  differences 
of  drivers  and  pedestrians  in  various  countries.  The 
element  of  individual  and  collective  behaviour  as 
measured  by  traffic  platoons  (dispersions,  headways  and 
spacing,  headway  stability)  either  constitutes  or 
significantly  influences  traffic  behaviour.  Traffic 
platoons,  therefore  is  one  measure  to  report  deviations 
from  the  accepted  level. 


11.  Preferential  street  treatment .  To  provide  guidelines  in 
the  traditional  preferential  street  treatment  of 
co-ordination  design.  Streets  with  compact  platoon 


. 
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movements  can  be  favoured  over  streets  with  considerable 
platoon  dispersion. 


12.  Ranking  of  network  links.  To  aid  in  the  ranking  of 
network  links  in  order  of  their  importance  in  terms  of 
effectiveness  of  co-ordination.  Link  ranking  is 
important  in  off-line  signal  control  strategies  such  as 
SIGRID  (SI_gnal  GRid  Design)  Program  (37)  where  links  are 
weighted  by  their  importance  before  optimization. 


13.  Need  for  coordination.  To  provide  an  indication  of 
co-ordination  needs.  For  example,  links  without  platoon 
movements  can  be  excluded  from  the  network  description 
to  save  computing  time. 


14.  Network  condensation.  To  aid  in  network  condensation  for 
easier  analysis.  For  example,  the  Combination 
Program(23)  which  is  only  applicable  to  certain  complex 
networks  can  be  condensed  to  suit  the  program  by  the 
deletion  of  links  with  relatively  ill-defined  platoon 
structures . 


15.  Understand i ng  the  flow  characteristics  at  a  midblock 
access  road  between  two  signal i zed  intersections. 

Estimating  the  shape  of  a  cyclic  flow  profile  flow  for  a 
two-way  road  at  a  midblock  access  capacity  analysis.  If 
the  platoon  from  both  directions  of  flow  cross  the 
midblock  access  road  at  the  same  time,  adequate  gaps  for 
the  minor  street  flow  may  be  created.  On  the  other  hand, 
staggering  the  arrival  of  both  platoons  may  imply  that  a 
capacity  deficiency  may  materialize  at  the  access  road. 


6 .  SUMMARY 


6 . 1  Conclusion 

The  results  of  this  research  have  shown  that  the 
dispersion  of  traffic  platoons  between  fixed-time  signalized 
intersections  can  be  simulated  quite  accurately  using 
Robertson's  Recurrence  Model.  In  the  absense  of  locally 
available  information,  Robertson's  recommended  empirical 

parameters  of  alpha=0.50  and  beta=0.80  or  more  recently 
alpha=0.44  and  beta=0.80  represent  good  starting 

approx i mat i ons .  Under  circumstances  such  as  winter  weather 
conditions  and  degrees  of  saturation  less  than  90%  on  the 
street  system,  Robertson's  recommended  parameters  will  be  in 
error.  The  effect  of  these  errors  has  a  signficant  effect  on 
the  accuracy  of  signal  offset  design  which  is  measured  by 
such  indicators  as  delays,  stops  and  performance  index. 
Significant  improvements  can  be  made  by  calibrating  the  two 
individual  empirical  parameters  and  introducing  an 

initialization  equation  (Modified  Robertson's  Recurrence 

Model)  to  start  the  prediction  process.  Specifically,  the 
following  errors  in  signal  performance  are  noted  : 

MINIMUM  VEHICULAR  DELAY 

1.  A  Robertson's  Recurrence  Model  which  has  not  been 
calibrated  to  local  conditions  could  be  in  error  up  to 
20%. 

2.  A  Robertson's  Recurrence  Model  which  has  been  calibrated 
to  local  conditions  could  be  in  error  up  to  16%. 
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3.  A  Modified  Robertson’s  Recurrence  Model  which  has  been 
calibrated  to  local  conditions  could  still  be  in  error 
up  to  7%. 


MINIMUM  VEHICULAR  STOPS 


1.  A  Robertson's  Recurrence  Model  which  has  not  been 
calibrated  to  local  conditions  could  be  in  error  up  to 
78%. 


2.  A  Robertson's  Recurrence  Model  which  has  been  calibrated 
to  local  conditions  could  be  in  error  up  to  43%. 


3.  A  Modified  Robertson's  Recurrence  Model  which  has  not 
been  calibrated  to  local  conditions  could  still  be  in 
error  up  to  38%. 


MINIMUM  PERFORMANCE  INDEX 


1.  A  Robertson's  Recurrence  Model  which  has  not  been 
calibrated  to  local  conditions  could  be  in  error  up  to 
14%. 


2.  A  Robertson's  Recurrence  Model  which  has  been  calibrated 
to  local  conditions  could  be  in  error  by  15%. 


3.  A  Modified  Robertson's  Recurrence  Model  which  has  been 
calibrated  to  local  conditions  could  still  be  in  error 
up  to  4%. 


Based 

on  the 

analy s i s 

of  six  survey 

sites,  the 

Researcher 

concludes 

that  both 

calibration  and 

modi f icat ions 

to  Robertson's  Recur 

rence  Model 

can  reduce 

the  error  in 

prediction  of  delay  by  a  factor  of  3,  error  in  prediction  of 
stops  by  a  factor  of  2,  and  error  in  prediction  of  the 
performance  index  by  a  factor  of  3.5.  These  conclusions  have 
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significant  implications  in  the  following  areas 


1  . 


The  design  of  signal  timings  along  a  one  wa 
where  the  objective  is  to  minimize  the  number 
An  error  in  the  estimate  of  the  number  of 
factor  of  2  will  result  in  sub-optimal  signal 


y  arterial 
of  stops, 
stops  by  a 
design . 


2.  The  design  of  network  timings  in  the  TRANSYT  model  is 
based  on  an  optimization  procedure.  The  offsets  between 
signals  are  adjusted  by  a  "hill-climbing"  optimization 
procedure  to  minimize  an  index  of  performance.  The  index 
of  performance  is  a  weighted  combination  of  delay  and 
stops.  Errors  in  delay  estimates  and  stops  estimates, 
when  combined  into  a  performance  index  generally  will 
result  on  average  in  a  larger  error.  Optimization  of  an 
index  which  is  in  error  will  result  in  sub-optimal 
network  timing  design. 


6.2  Recommendations 

To  improve  on  the  prediction  of  several  signal 
performance  indicators  such  as  delay,  stops  and  performance 
index  the  following  recommendations  are  put  forth: 

1.  Robertson's  Recurrence  Model  in  the  TRANSYT  method  of 
signal  design  should  be  modified  by  introducing  an 
initialization  equation  to  the  prediction  process. 


2.  Basic  imput  into  the  TRANSYT  method  of  signal  design 
should  be  modified  to  allow  the  individual  empirical 
parameters  of  alpha  and  beta  as  input  parameters. 

3.  In  the  absence  of  data  for  a  specific  signal  network, 
the  following  empirical  parameters  should  be  used  in 
Edmonton,  Alberta  (Canada): 


SUMMER  CONDITIONS  :  Alpha  =  0.34,  Beta  =  0.75 
WINTER  CONDITIONS  :  Alpha  =  0.34,  Beta  =  0.66 


4.  The  Micro-computer  Cyclic  Flow  Profile  Monitor  has  been 
demonstrated  to  be  an  excellent  tool  for  collecting 
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real-time  traffic  platoon  data.  It  is  estimated  that 
with  two  technical  support  staff  and  two-three  hours  of 
total  survey  time  (setting  up  of  the  machine,  survey 
itself,  dis-assembly  of  survey  equipment),  approximately 
one-half  to  one  hour  of  data  can  be  collected  to  a 
resolution  of  one-tenth  of  a  second.  With  the  software 
packages  ( BRAD#+*PLOTLID ,  GORDONV 1 ,  GORDONV2  and  DSPI ) 
that  were  specifically  designed  for  this  research,  it  is 
estimated  that  one  to  one-half  hour  of  support  staff 
time  would  be  required  to  obtain  the  most  appropriate 
value  of  alpha  or  beta.  The  City  Of  Edmonton  traffic 
signal  design/analysis  team  is,  therefore,  encouraged  to 
use  the  data  collection  technique  and  the  various 
software  programs  developed  for  this  research.  This  on 
going  process  of  data  collection  and  analysis  at 
selective  roadway  links  i s  needed  to  confirm  that  the 
predict  ions  by  TRANSYT  reflect  the  actual  field 
conditions . 


6.3  Further  Research  Needs 

The  following  further  research  needs  have  been 
ident i f ied : 


The  initialization  equation  introduced  into  the  Modified 
Recurrence  Model  was  based  on  the  assumption  of 
steady-state  conditions.  Additional  research  into  the 
validity  of  this  assumption  is  required  especially  when 
it  is  observed  during  peak  time  periods,  that  time 
varying  demands  occur. 


2.  The  analysis  of  the  data  base  presented,  indicates  that 
the  individual  value  of  the  alpha  and  beta  parameters  is 
influenced  by  weather  conditions  (summer  vs  winter)  and 
degree  of  saturation.  Specifically,  it  was  observed 
that,  with  the  summer  data,  a  consistent  trend  (as 
degree  of  saturation  increased  the  individual  value  of 
alpha  and  beta  increased)  appeared.  More  winter  and 
summer  data  could  be  instrumental  in  defining  more 
refined  design/  analysis  guidelines  in  the  use  of  the 
two  empirical  parameters. 


3.  More  study  is  needed  to  explain  why  a  Calibrated 
Modified  Robertson’s  Recurrence  Model  still  results  in 
an  underestimation  of  minimum  delay  by  7%,  an 
underestimation  of  minimum  stops  by  38%,  and  an 
underestimation  of  minimum  performance  index  by  4%. 
Suggestions  include  :  revisions  to  the  form  of  the  F 
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factor  (making  /?T  a  real  value  vs  an  integer  value), 
prediction  downstream  is  based  on  a  portion  of  the 
previous  calculated  interval  (i+^T-1  vs  i+/JT-2  (  3 , 4  .  .  .  )  )  . 


4.  Detail  statistical  analysis  of  the  travel  time  of 
individual  vehicles  in  a  platoon  (platoon  followers)  vs 
the  travel  time  of  platoon  leaders  so  that  a  better 
understanding  of  the  beta  empirical  parameter  can  be 
determined . 
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